Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,145 G→C 6.2% intergenic (+157/+33) sibC → / ← serA sRNA antisense regulator of toxic IbsC protein/D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,057,1450GC 6.2% 245.3 / 7.0 92intergenic (+157/+33)sibC/serAsRNA antisense regulator of toxic IbsC protein/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (33/53);  new base C (3/3);  total (36/56)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGG  >  NC_000913/3056998‑3057289
                                                                                                                                                    |                                                                                                                                                
gCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGt                                                                                                                                                 <  2:837169/149‑1 (MQ=255)
  cccATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGa                                                                                                                                               >  2:55316/1‑149 (MQ=255)
   ccATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAATGGATAAAAAAACGG‑GCAAGTCAGTGAc                                                                                                                                              <  1:544769/149‑1 (MQ=255)
   ccATGCTGCCCTGCCATCGTAAAACCCCATAAAATAACCACAACGCATTGATCTGACTTTGATTTATTTTCTGCAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGAc                                                                                                                                              <  2:904336/149‑1 (MQ=255)
         tgccctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGCATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGccc                                                                                                                                        >  2:762463/1‑149 (MQ=255)
         tgccctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAGGTCACTGACTTGccc                                                                                                                                        <  1:345458/149‑1 (MQ=255)
         tgccctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGccc                                                                                                                                        <  2:1012172/149‑1 (MQ=255)
           ccctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAAAGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGt                                                                                                                                      >  1:174174/1‑149 (MQ=255)
           ccctgccATCGAAAATCCCCATTAAATAAACACCACGCATTGATATGACATTGATTTATATTCTGCAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGt                                                                                                                                      <  1:774610/149‑1 (MQ=255)
               gccATCGTAAATCCCCATTAAATAAACACACCGCCTTGCTCTGACTTTGCTTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGGGGAGCAGGGATAAAAAAACGG‑GCAGGGCACTGACTTGCCCGttttt                                                                                                                                  <  1:755397/149‑4 (MQ=255)
               gccATCGTAAATCCCCATTAAATAAACACAACGCAGTGATCTGACTTTGAGTTATTTTCTGGAGCAGACTCGCAAAGGAGAATGCGCACCGCGGCAACGGTGAGGAGAAGGGATAAAAAAACGG‑GCAGGGCACTGACTTGCCCGttttt                                                                                                                                  <  1:498620/149‑4 (MQ=255)
                 cATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGCCCTGCCCGTTGAttt                                                                                                                                >  2:115788/1‑149 (MQ=255)
                 cATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGAttt                                                                                                                                >  2:153022/1‑149 (MQ=255)
                  aTCGTAAATCCCCATTAAATAATCACAACGCATTGATCTGACATTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCATCGCGGCAACTGTGTGGAGAACTGTTAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGAttgt                                                                                                                               >  1:697048/1‑147 (MQ=255)
                  aTCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAg                                                                                                                                                  >  2:375982/1‑130 (MQ=255)
                  aTCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAg                                                                                                                                                  <  1:375982/130‑1 (MQ=255)
                    cGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCa                                                                                                                             <  2:380281/149‑1 (MQ=255)
                    cGTAAATCCCCATTAAATAAACACAACGCATTGATCTGAAGTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCa                                                                                                                             <  2:582180/149‑1 (MQ=255)
                     ttAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                                                                            <  1:1173850/148‑1 (MQ=255)
                     gTAAATCCCCATTAAAAAACCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGACTGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                                                                            <  1:1178262/149‑1 (MQ=255)
                           ccccATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                                                                      <  2:48575/149‑1 (MQ=255)
                            cccATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAggg                                                                                                                     <  1:1099977/149‑1 (MQ=255)
                             ccATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCAGTTGATTTTCAGAGAAgggg                                                                                                                    >  2:14375/1‑149 (MQ=255)
                             ccATTAAATAAACACAACGCATTGATCCGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                                                                    <  1:1132006/149‑1 (MQ=255)
                              cATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGa                                                                                                                   <  2:387339/149‑1 (MQ=255)
                                    ataaaCACAACGCATTTATCTGACTTTGATTTATTTTCTGGAGCAGACACGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGCAAAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                                                             >  1:552806/1‑149 (MQ=255)
                                       aaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAca                                                                                                          <  2:1034508/149‑1 (MQ=255)
                                        acacAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAccg                                                                                                         >  2:220731/1‑147 (MQ=255)
                                        acacAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                                                                         >  1:1035691/1‑149 (MQ=255)
                                         aacaACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGTATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagc                                                                                                        <  1:816529/148‑1 (MQ=255)
                                                cATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACggg                                                                                                 <  1:228462/149‑1 (MQ=255)
                                                 aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgc                                                                                                >  1:388220/1‑149 (MQ=255)
                                                   tGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcgc                                                                                              <  1:390689/149‑1 (MQ=255)
                                                      tCTGACTTTGATTTTTTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACTG‑TCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGTGGAATTAGTACAGCAGACGGGCGCGaa                                                                                           <  2:495794/149‑1 (MQ=255)
                                                          aCTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGt                                                                                       >  1:1122259/1‑149 (MQ=255)
                                                            tttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAc                                                                                                           >  1:565853/1‑127 (MQ=255)
                                                            tttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAc                                                                                                           <  2:565853/127‑1 (MQ=255)
                                                            tttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAc                                                                                     >  1:1150512/1‑149 (MQ=255)
                                                             tttatttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAcc                                                                                    <  1:471786/146‑1 (MQ=255)
                                                              tGATGTATTTTCTCGAGCAGACTCCCAAAGTAGAATCCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAATCGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGCGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAccc                                                                                   <  2:54461/149‑1 (MQ=255)
                                                               tatttatttTCTGGAGCAGACTCGCAAAGTCGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGTAATTAGTACAGCAGACGGGCGCGAATGGTACCCg                                                                                  <  2:68200/148‑1 (MQ=255)
                                                               gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCg                                                                                  <  2:1122259/149‑1 (MQ=255)
                                                               gATTTATTTTCTGGAGCAGACTCGCAAAGAAGAATGCGCAACGCGGCAACGGGGTGGAGAAGGGATAAAAAAACGG‑GCAAGGCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCg                                                                                  <  1:288338/149‑1 (MQ=255)
                                                                atttatttTCTGGAGCAGCCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGTTAAAAAAACGG‑GCAAGTCAGTGACCTGCGCGGTGATTTTCAGAGAAGGGGATGTAGTACAGCAGCCGGGCGCGAATGGTACCCgg                                                                                 >  1:880129/1‑149 (MQ=255)
                                                                atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAATCGG‑GCAAGTCAGTGCCCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGGAGACGGGCGCGAATGGTACCCgg                                                                                 >  1:1061903/1‑149 (MQ=255)
                                                                atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCgg                                                                                 >  2:236018/1‑149 (MQ=255)
                                                                            gAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAt                                                                     >  1:125516/1‑149 (MQ=255)
                                                                             atcAGACTCGCAAAGTAGACTGCTCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGAACTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAtt                                                                    <  2:813921/147‑1 (MQ=255)
                                                                               cAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                                                                  <  2:692557/149‑1 (MQ=255)
                                                                               cAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCAGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                                                                  >  2:760018/1‑149 (MQ=255)
                                                                               cAGACGCTCAACGTAGAATGCGCAACGCGGCAACGGTGTGGATAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                                                                  <  2:293957/149‑1 (MQ=255)
                                                                                aGACTCGCATAGTAGAATGCGCACCGCGGCATCGGTGTGGAGATGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGTCGCGAATGGTACCCGGTATAGCTTTCATTGcc                                                                 <  1:543868/149‑1 (MQ=255)
                                                                                aGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGAGTGGAGAAGTGATAAAAAAACGC‑TCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc                                                                 <  1:672406/149‑1 (MQ=255)
                                                                                 gACTCGCAAAGTAGCATGCGCAACGCGGCAACGGTGTTGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATATTCAGAGATGGGGAATTAGTACAGGAGACGGGCGCGTATGGTAACCGGAATAGCTTTCATTGCCt                                                                >  1:661972/1‑149 (MQ=255)
                                                                                 caCTCGCACAGGAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAACCGG‑GCACGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt                                                                <  1:852388/148‑1 (MQ=255)
                                                                                  aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGATGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGTTTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTg                                                               <  1:770677/149‑1 (MQ=255)
                                                                                  aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTg                                                               <  1:55316/149‑1 (MQ=255)
                                                                                  aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTg                                                               <  1:492924/149‑1 (MQ=255)
                                                                                   ctctCAAAGTAGAATGCGCATCGCGGCAACGTTGTGGAGTAGGGATTAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGTAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc                                                              <  1:1186117/149‑1 (MQ=255)
                                                                                    tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                                                             <  2:469819/149‑1 (MQ=255)
                                                                                      gCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc                                                           <  2:1035691/149‑1 (MQ=255)
                                                                                         aaGTAGAATGCGCAACGCGGCAACCGTGTGGACAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt                                                        <  2:667678/149‑1 (MQ=255)
                                                                                              gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCCGCAGCGCTTTTTc                                                   <  1:414728/149‑1 (MQ=255)
                                                                                               aaggtgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAACCGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCg                                                  <  1:521460/144‑1 (MQ=255)
                                                                                               aaTGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGATTTAGTACAGCAGACGGGCGCGAATGGAACCCGGAATAGCTTTCATTGCCTGCAGCGCTTGTTCg                                                  >  2:567772/1‑149 (MQ=255)
                                                                                               aaTGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCg                                                  >  2:563322/1‑149 (MQ=255)
                                                                                                   cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCaa                                              >  1:932923/1‑149 (MQ=255)
                                                                                                    gcAACGCGGCAACGGGGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTGCAGAGAAGGTGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGTCTTAGCTTTCTTTGCCTGCTGCGCTTTTTCGGCAAc                                             >  1:666176/1‑149 (MQ=255)
                                                                                                     cAACGCGGCAACTGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACg                                            <  2:748500/149‑1 (MQ=255)
                                                                                                        cgcgGCAAGGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCt                                         <  1:335637/149‑1 (MQ=255)
                                                                                                           ggCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTTCCTGCAGCGCTTTTTCGGCAACGTCTtcg                                      >  1:943361/1‑149 (MQ=255)
                                                                                                            gCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTACAGAGAAGGGGAATTAGTACAGCAGACGGGCCCGAATGGAACCCGGAATAGCTTTCATTGCCTGCAGCGATTATTCGGCAACGACTtcgt                                     >  2:563561/1‑149 (MQ=255)
                                                                                                                  gtgtGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGTAATTTGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCtt                               <  2:174174/149‑1 (MQ=255)
                                                                                                                   tgtgGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCTGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCtta                              >  2:1006504/1‑148 (MQ=255)
                                                                                                                      ggAGAACGGATAAAAAAACGG‑GCAAGTCAGTGACGTGGCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGGGCTTTTTCGGCAACGTCTTCGTCGGCTtcaat                           <  1:800325/149‑1 (MQ=255)
                                                                                                                       gagaAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaata                          <  1:742683/149‑1 (MQ=255)
                                                                                                                         gaAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatc                        <  1:153022/149‑1 (MQ=255)
                                                                                                                           aGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaa                      <  1:14375/149‑1 (MQ=255)
                                                                                                                            gggATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaat                     >  2:52088/1‑149 (MQ=255)
                                                                                                                            gggATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaat                     <  2:552806/149‑1 (MQ=255)
                                                                                                                            gggATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTGCGTCGGCTtcaatatcaat                     >  2:1098501/1‑149 (MQ=255)
                                                                                                                                aaaaaaaaCGGCGCAAGGCCGTCACCTGCCCGTTGATTTTCAGAGCAGGGGAATTAGAACAGCAGACGGGCGCGATTGGTTCCCGGAATAGCTTTCATTGGCTGCAGCGCTTTTTCGGCATCGGCTTCGTCGGCTTCAATATCAATAAc                  <  2:333038/148‑1 (MQ=255)
                                                                                                                                   aaaaaCGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACcaca              >  2:866704/1‑149 (MQ=255)
                                                                                                                                     caaCGG‑GCAGGTCACTGACTTGCACGTTGATTTTCAGAGAAGGGGAATAACTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCACCGTATTCGTCGGCTTCAATATCAATAACcacaga            >  1:241352/2‑147 (MQ=255)
                                                                                                                                      aaCGG‑GCAAGTCAGAGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGCCGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCGTCAATATCAATAACCACATaa           >  1:55943/1‑149 (MQ=255)
                                                                                                                                       aCGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAAc          <  2:152446/149‑1 (MQ=255)
                                                                                                                                         gg‑gCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAAccc        <  2:2617/149‑1 (MQ=255)
                                                                                                                                            gCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCAt      >  2:1109839/1‑149 (MQ=255)
                                                                                                                                            gCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCGTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATTACCCAt      >  2:1105855/1‑149 (MQ=255)
                                                                                                                                             cAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATc     >  1:1035508/1‑149 (MQ=255)
                                                                                                                                               aGTCAGTGACCGGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTg   <  2:1087845/149‑1 (MQ=255)
                                                                                                                                               aGTCAGTGACCAGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTg   <  1:760018/149‑1 (MQ=255)
                                                                                                                                                gTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTgg  <  1:931419/149‑1 (MQ=255)
                                                                                                                                                gTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCATCAATATCAATAACCACATAACCCATCTgg  >  2:195517/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGG‑GCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGG  >  NC_000913/3056998‑3057289

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑