Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,193,286 G→C 5.8% D97E (GAC→GAG thiF ← adenylyltransferase, modifies ThiS C‑terminus

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,2860GC 5.8% 359.2 / 9.7 127D97E (GAC→GAGthiFadenylyltransferase, modifies ThiS C‑terminus
Reads supporting (aligned to +/- strand):  ref base G (48/71);  new base C (4/4);  total (52/75)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCG  >  NC_000913/4193141‑4193434
                                                                                                                                                 |                                                                                                                                                    
gcggcgTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCg                                                                                                                                                   <  1:1412668/149‑1 (MQ=255)
  ggcgTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGCGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCggg                                                                                                                                                 <  2:283931/149‑1 (MQ=255)
   gcgTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGa                                                                                                                                                >  2:1337537/1‑149 (MQ=255)
        ttGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCaa                                                                                                                                           <  2:211670/149‑1 (MQ=255)
          gAGTGCCACGCATGAGGCTTTTATCTCGTGGGGAGTCGCCATATTGTCGGTACATTCGAGCACCAGATCGGCCCGTGCAACCGCATCTTTTTACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTGAGTTGAATGTCGGGATTCAACt                                                                                                                                         <  2:1462081/149‑1 (MQ=255)
          gAGTGCCACGCAGGCGGCATTAATCTCCTGGCGCGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATCCCGACATTCAACt                                                                                                                                         <  2:1053886/149‑1 (MQ=255)
          gAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGGCGGGATTCAACt                                                                                                                                         <  1:1208253/149‑1 (MQ=255)
          gAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATCCCGACATTCAACt                                                                                                                                         <  1:1571357/149‑1 (MQ=255)
          gAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTAGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATCCCGACATTCAACt                                                                                                                                         <  2:737507/149‑1 (MQ=255)
            gTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACtgt                                                                                                                                       <  1:1233958/149‑1 (MQ=255)
            gTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACtgt                                                                                                                                       <  1:395624/149‑1 (MQ=255)
              gCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTgtgt                                                                                                                                     <  2:90828/149‑1 (MQ=255)
              gCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTgtgt                                                                                                                                     <  1:848447/149‑1 (MQ=255)
                   gCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTc                                                                                                                                <  2:472681/149‑1 (MQ=255)
                     aGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCTCATCGGCCCGTGCAACCGCATCTTTTAACGCCTCTCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGc                                                                                                                              >  1:139802/1‑149 (MQ=255)
                     aGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTACCGCCTCCCCCGTTAACAGTTGTTGTAATGCTGTCAGTTGAATCGCGGGAGGCATCTGTGTCAGTCGc                                                                                                                              >  2:851161/1‑149 (MQ=255)
                         ggcATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGtt                                                                                                                          <  2:204795/149‑1 (MQ=255)
                         ggcATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGtt                                                                                                                          <  1:795867/149‑1 (MQ=255)
                          gcATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACTGGCTAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTATCCGTTGTTGTTTTGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCTGTCGCTGTTg                                                                                                                         <  1:925372/149‑1 (MQ=255)
                          gcATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTg                                                                                                                         >  1:759676/1‑149 (MQ=255)
                           cATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCCCCCGTTAACGGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCGGTTgg                                                                                                                        >  1:791158/1‑149 (MQ=255)
                           cATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTgg                                                                                                                        <  1:1401863/149‑1 (MQ=255)
                             ttAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCt                                                                                                                      <  1:1386112/149‑1 (MQ=255)
                                aTCTCCTGGCGAGTCGCCATATTGTCGGTTCAGTCGATCTCCACATCGGCCCGTGCAACCGCATCTTTTTACGCCTCACGCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTAGGGATTCAACTGTGTCAGTCGCTGTTGGCTGAc                                                                                                                   <  2:1080071/149‑1 (MQ=255)
                                aTCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCCCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGAc                                                                                                                   >  1:236382/1‑149 (MQ=255)
                                    ccTGGCGAGTCGCCATTTTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACAGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGc                                                                                                               <  2:1002472/149‑1 (MQ=255)
                                    ccTGGCGAGGCGCCATATTGTAGGTACAGTCGAGCACCACCTCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTACCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGc                                                                                                               <  2:911763/149‑1 (MQ=255)
                                      tGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGACTCCCCCGTTAACCGTTCTTGTAATGCTGTCAGTTGAATGTCGGGATACAAATGTGTCAGTCGCTGTTGGCTGCCCTGCGa                                                                                                             >  2:366662/1‑149 (MQ=255)
                                        gCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACAACGGCCCGTGCAAACGCATCTATTAACGCCTCACCCGTTAACCCATGTAGTAATGCTTTCAGTTTAATGTCTGGATACAACTGTGTCAGTAGCTGTTGGCTGACCGGCGAtt                                                                                                           >  2:1184207/1‑149 (MQ=255)
                                           aGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGGTGTAATGCTGCCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCg                                                                                                        <  1:548505/149‑1 (MQ=255)
                                            gTCGCCATATTGTCGGGACAGTCGATGACCACATCGGCCCGTGCAACCGCATCTTTTATCGCCTCACCCGTTAACCTTTGTTGTAATGCTGTCTGTTGAATGTCGGGATTCATCTGTGTCAGTCGCTGTTTGCTGACCTGCGATTTCgg                                                                                                       <  2:1427752/149‑1 (MQ=255)
                                             tCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCgg                                                                                                       >  1:309861/1‑148 (MQ=255)
                                             tCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCgg                                                                                                       <  2:309861/148‑1 (MQ=255)
                                              cGCCATATTGTATGTACAGTCGAGCACCACATCGTCCCGTGCATCCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGc                                                                                                     <  1:980335/149‑1 (MQ=255)
                                                 cATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgat                                                                                                  <  2:1006052/149‑1 (MQ=255)
                                                 cATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgat                                                                                                  >  1:1502123/1‑149 (MQ=255)
                                                    atTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTTTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatcga                                                                                               <  1:1563759/149‑1 (MQ=255)
                                                     ttGTCGGTACAGTCGCGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatcgat                                                                                              <  1:1062469/149‑1 (MQ=255)
                                                      tGTCGGTACATTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGata                                                                                             <  1:911492/149‑1 (MQ=255)
                                                      tGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCTCCCGGTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGATGACCTGCGATATCGGGCGATCGata                                                                                             <  1:1496849/149‑1 (MQ=255)
                                                        tCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATc                                                                                           >  1:56924/1‑149 (MQ=255)
                                                        tCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATc                                                                                           <  1:211207/149‑1 (MQ=255)
                                                           gTACAGTCGAGCACCACATCGGCCCGTGCACCCGCCTCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGCCCTGCGATTTCGGGCGATCGATttcttc                                                                                        >  2:1185507/1‑149 (MQ=255)
                                                           gTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTc                                                                                        >  2:427393/1‑149 (MQ=255)
                                                            tACAGTCGAGCACCACATCGGCCCGTGCACCCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCa                                                                                       >  1:220038/1‑149 (MQ=255)
                                                               aGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCGTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCACCAGTGTCCGTCGCTGTTGGCTGCCCTGCTATTTCGGGCGATCGATATCTTCAGTg                                                                                    >  1:288122/1‑149 (MQ=255)
                                                                gTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTgg                                                                                   >  1:542509/1‑149 (MQ=255)
                                                                gTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTgg                                                                                   >  1:527586/1‑149 (MQ=255)
                                                                    aGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTaaa                                                                               >  2:1102353/1‑149 (MQ=255)
                                                                      caccacATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAaga                                                                             <  2:1359628/149‑1 (MQ=255)
                                                                        ccacATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAgg                                                                           >  2:561949/1‑149 (MQ=255)
                                                                        ccacATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAgg                                                                           <  1:985818/149‑1 (MQ=255)
                                                                        ccacATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAgg                                                                           >  1:1098707/1‑149 (MQ=255)
                                                                        ccacATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAgg                                                                           <  1:280950/149‑1 (MQ=255)
                                                                        ccacATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGGCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGAGATCTTCAGTGGTAAAGAgg                                                                           >  2:153187/1‑149 (MQ=255)
                                                                        ccacATCGGCCCGGGCAACCGCATCTTTTAACGCCTCACCCTTTAACCGATGTTGTAATGCTGTCAGTTGAATGTCGGTATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGAGGAAAAGAgg                                                                           <  2:1048961/149‑1 (MQ=255)
                                                                         cacaTCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGCTATCTTCAGTGGTAAAGAGGa                                                                          >  2:57556/1‑149 (MQ=255)
                                                                         cacaTCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGa                                                                          <  1:513805/149‑1 (MQ=255)
                                                                           caTCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGAtt                                                                        >  2:224950/1‑149 (MQ=255)
                                                                            aTCGGCCCGTGCAACCGCAGCTATTACCGCCTCACCCGTTAACCGTTGTTGTAATGCTGGAAGTAGAATGTCGGGAGTCAACTGTGCCAGGCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGAttt                                                                       <  2:339025/149‑1 (MQ=255)
                                                                              cGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGt                                                                     <  2:1124333/149‑1 (MQ=255)
                                                                                gCCCGTGCAACCGCATCTTTTACCGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCg                                                                   >  1:1051377/1‑149 (MQ=255)
                                                                                 cccGTGCAACCGCATCTTTTAACGCCTCACCCGTTACCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGt                                                                  <  1:148309/149‑1 (MQ=255)
                                                                                 cccGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGt                                                                  <  2:423648/149‑1 (MQ=255)
                                                                                   cGTGCAACCGCATCTTTTAACGCCTCACCCGTAAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTg                                                                <  1:1192565/149‑1 (MQ=255)
                                                                                     tGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCa                                                              <  1:902346/149‑1 (MQ=255)
                                                                                        aaCCGCATCTTTTAACGCCTCACCCGTTAACCTTTGTTGTAATGCTGTCAGTTCAATGTCTGGATTCAACTCTGTCAGTCGCTGATGGCTGACCTGCGATTTCGGGCGATCGATATCGTCAGTGGTAAAGAGTATTTGTCGTTGCACAt                                                           >  1:1527337/1‑149 (MQ=255)
                                                                                        aaCCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCCGTCGCGGTTGGCTGACCGGCGATTTCGGGCGATCTATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGAt                                                           >  2:371889/1‑149 (MQ=255)
                                                                                        aaCCGCATCTTTTAACGCCTCACCCGTTAACCGGTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCCACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGAGTTGTCGTTGCAGAt                                                           >  2:779021/1‑149 (MQ=255)
                                                                                         aCCGCATCTTTTAACGCCTCACCCGTTATCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGAtt                                                          <  1:1253123/149‑1 (MQ=255)
                                                                                          ccGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTg                                                         >  1:475575/1‑149 (MQ=255)
                                                                                          ccGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTg                                                         <  1:903748/149‑1 (MQ=255)
                                                                                          ccGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGCAAAGAGGATTTGTCGTTGCAGATTg                                                         <  1:401589/149‑1 (MQ=255)
                                                                                            gCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCt                                                       <  2:1230029/149‑1 (MQ=255)
                                                                                             cATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCtt                                                      >  1:734629/1‑149 (MQ=255)
                                                                                               tCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGGCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTaa                                                    <  1:366662/149‑1 (MQ=255)
                                                                                                 ttttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAAt                                                  <  1:538623/149‑1 (MQ=255)
                                                                                                  tttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATg                                                 >  1:1240436/1‑149 (MQ=255)
                                                                                                  tttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATg                                                 <  1:1279957/149‑1 (MQ=255)
                                                                                                  tttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATg                                                 >  2:908305/1‑149 (MQ=255)
                                                                                                  tttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCACTGGTAAAGAGGATTTGTGGTTGCAGATTGCTTAAATg                                                 >  1:851260/1‑149 (MQ=255)
                                                                                                  tttAACGCCTCACCCGGTAACCGGTGTTGTCATGCTGTCAGTTGACTGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATg                                                 <  2:328125/149‑1 (MQ=255)
                                                                                                     aaCGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGcac                                              <  2:225795/149‑1 (MQ=255)
                                                                                                         ccTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGCTCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcg                                          >  2:1026753/1‑149 (MQ=255)
                                                                                                         ccTCACCCGTAAACCGTTGTTGTAATGCTGTCCGGTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcg                                          <  2:791158/149‑1 (MQ=255)
                                                                                                          cTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGCCCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATTCAAAtctt                                         >  2:1133516/1‑147 (MQ=255)
                                                                                                           tCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGTCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcgta                                        >  2:650495/1‑148 (MQ=255)
                                                                                                             aCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATcgtcgt                                      >  2:240409/1‑149 (MQ=255)
                                                                                                                 gTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGAGTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATc                                  <  1:1541488/149‑1 (MQ=255)
                                                                                                                    aaCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGTATTCAACTGTGTCAGGCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGc                               <  1:224950/149‑1 (MQ=255)
                                                                                                                    aaCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGTCGATCGATATCTTCAGTGGTAAAGAGGATTTGGCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCAGc                               <  1:842774/149‑1 (MQ=255)
                                                                                                                      ccGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCa                             <  1:513667/149‑1 (MQ=255)
                                                                                                                       cGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAg                            <  1:600733/149‑1 (MQ=255)
                                                                                                                            ttgtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCACTCGCTGTTTGCTGACCTGCGATTTCGGGCGATCGATATCTACAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCAAAACTTCGTCATCTGACAGTAACa                       >  1:349576/1‑149 (MQ=255)
                                                                                                                              gtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGc                     >  1:1314231/1‑149 (MQ=255)
                                                                                                                                      ttCAGTTGTATGTCGGGATTCAACTGTGTGAGTCGCTGTTGTCTGACCTGCTATTTCGGGCGATGGAAATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGGCCCGAc             <  2:1471323/148‑1 (MQ=255)
                                                                                                                                       acAGTTGAATCCCGACATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACg            >  2:1327936/2‑149 (MQ=255)
                                                                                                                                            tGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATATCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGGCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAgc       <  1:1083890/149‑1 (MQ=255)
                                                                                                                                                tGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGcgccc   >  1:624704/1‑149 (MQ=255)
                                                                                                                                                 gTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGAcgcccgcgcccg  >  2:392777/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCG  >  NC_000913/4193141‑4193434

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: