Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,125,308 C→T 5.5% R51H (CGC→CAC)  glcF ← glycolate dehydrogenase, putative iron‑sulfur subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,125,3080CT 5.5% 226.2 / 2.1 72R51H (CGC→CAC) glcFglycolate dehydrogenase, putative iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base C (29/39);  new base T (2/2);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGG  >  NC_000913/3125160‑3125452
                                                                                                                                                    |                                                                                                                                                
atCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATgcg                                                                                                                                                 >  1:462502/1‑150 (MQ=255)
atCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATgcg                                                                                                                                                 >  1:497713/1‑150 (MQ=255)
atCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATgcg                                                                                                                                                 >  1:700021/1‑150 (MQ=255)
  cACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGTGcc                                                                                                                                               >  1:745843/1‑150 (MQ=255)
       ccGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCcgcgc                                                                                                                                          <  1:306899/150‑1 (MQ=255)
          atatCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGc                                                                                                                                       <  2:549048/150‑1 (MQ=255)
           tatCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGcc                                                                                                                                      <  2:1479579/150‑1 (MQ=255)
               cAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTc                                                                                                                                  >  1:1508002/1‑150 (MQ=255)
                aGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTcc                                                                                                                                 >  2:1139570/1‑150 (MQ=255)
                 gCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCa                                                                                                                                <  1:1203704/150‑1 (MQ=255)
                            tAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGTGCCCGCGCGGCCCGTCCAGTTCATCGccc                                                                                                                     >  1:1408207/1‑150 (MQ=255)
                            tAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGccc                                                                                                                     >  2:1428695/1‑150 (MQ=255)
                            tAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGccc                                                                                                                     <  2:811180/150‑1 (MQ=255)
                             aGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCa                                                                                                                    <  2:1254056/150‑1 (MQ=255)
                                     ccAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGa                                                                                                            <  1:163241/150‑1 (MQ=255)
                                       aGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATa                                                                                                          >  1:95329/1‑150 (MQ=255)
                                       aGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATa                                                                                                          >  1:101558/1‑150 (MQ=255)
                                           ggACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGtt                                                                                                      <  2:1346902/150‑1 (MQ=255)
                                               aGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGc                                                                                                  <  1:853144/150‑1 (MQ=255)
                                                ggtggtTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCa                                                                                                 <  1:1402029/150‑1 (MQ=255)
                                                   ggtTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGt                                                                                              <  2:462502/150‑1 (MQ=255)
                                                        cacaATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCgg                                                                                         <  1:152869/150‑1 (MQ=255)
                                                        cacaATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCgg                                                                                         <  1:1428695/150‑1 (MQ=255)
                                                         acaATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTGCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGt                                                                                        <  1:1275584/150‑1 (MQ=255)
                                                              tACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACaaa                                                                                   <  1:301016/150‑1 (MQ=255)
                                                               aCGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACaaaa                                                                                  <  2:137978/150‑1 (MQ=255)
                                                               aCGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCAAGAAGCTGATAGGTTGGGCACGTTGCGGTACaaaa                                                                                  >  2:1548684/1‑150 (MQ=255)
                                                                cGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAAt                                                                                 >  1:906290/1‑150 (MQ=255)
                                                                 ggCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATc                                                                                <  2:1645926/150‑1 (MQ=255)
                                                                  gCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATcc                                                                               <  1:868610/150‑1 (MQ=255)
                                                                   cAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCg                                                                              <  1:1628215/150‑1 (MQ=255)
                                                                    aaGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGc                                                                             <  1:761420/150‑1 (MQ=255)
                                                                      gTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAg                                                                           <  1:210594/150‑1 (MQ=255)
                                                                        gAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTg                                                                         <  2:1659849/150‑1 (MQ=255)
                                                                        gAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTg                                                                         >  1:958890/1‑150 (MQ=255)
                                                                           gCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAac                                                                      >  2:1200676/1‑150 (MQ=255)
                                                                           gCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAac                                                                      >  1:1524607/1‑150 (MQ=255)
                                                                           gCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAac                                                                      <  1:875357/150‑1 (MQ=255)
                                                                            cAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAaca                                                                     >  2:1406292/1‑150 (MQ=255)
                                                                                gATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAgg                                                                 <  1:365556/150‑1 (MQ=255)
                                                                                gATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAgg                                                                 >  2:116266/1‑150 (MQ=255)
                                                                                gATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAgg                                                                 >  1:1532020/1‑150 (MQ=255)
                                                                                  tCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGgcg                                                               >  2:127659/1‑150 (MQ=255)
                                                                                      gaTGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGTGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCa                                                           <  2:745843/150‑1 (MQ=255)
                                                                                        tGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAgg                                                         >  2:112105/1‑150 (MQ=255)
                                                                                        tGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAgg                                                         >  1:283497/1‑150 (MQ=255)
                                                                                        tGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAgg                                                         <  2:1328642/150‑1 (MQ=255)
                                                                                                   tttAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGc                                              <  1:1616443/150‑1 (MQ=255)
                                                                                                   tttAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGc                                              >  1:317613/1‑150 (MQ=255)
                                                                                                   tttAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGc                                              <  2:1539732/150‑1 (MQ=255)
                                                                                                      aaGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTc                                           <  1:127659/150‑1 (MQ=255)
                                                                                                         cGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAg                                        >  1:778413/1‑150 (MQ=255)
                                                                                                           tGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAgcg                                      <  1:1406292/150‑1 (MQ=255)
                                                                                                               ttCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCTGGAAGCCGACAGCATCCTgcgcgcc                                  <  1:1188835/150‑2 (MQ=255)
                                                                                                               ttCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAgcgcgcg                                  >  2:281238/1‑150 (MQ=255)
                                                                                                                 cGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAgcgcgcgcg                                <  2:775373/150‑1 (MQ=255)
                                                                                                                    tGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGt                             >  1:631519/1‑150 (MQ=255)
                                                                                                                      ccTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTg                                                                         >  1:353379/1‑104 (MQ=255)
                                                                                                                      ccTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTg                                                                         <  2:353379/104‑1 (MQ=255)
                                                                                                                      ccTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTc                           <  2:1508002/150‑1 (MQ=255)
                                                                                                                         tCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGc                        <  2:497713/150‑1 (MQ=255)
                                                                                                                              cACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCAt                   >  2:90503/1‑150 (MQ=255)
                                                                                                                                ccTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCAtct                 >  1:32344/1‑150 (MQ=255)
                                                                                                                                 cTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCAtctc                >  1:1191133/1‑150 (MQ=255)
                                                                                                                                  tGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATctct               >  2:842924/1‑150 (MQ=255)
                                                                                                                                   gTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCtt              >  2:667831/1‑150 (MQ=255)
                                                                                                                                    tttAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTc             <  2:631519/150‑1 (MQ=255)
                                                                                                                                      tAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAg           >  2:1121985/1‑150 (MQ=255)
                                                                                                                                      tAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAg           >  1:366499/1‑150 (MQ=255)
                                                                                                                                       aaTCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGt          <  2:700021/150‑1 (MQ=255)
                                                                                                                                       aaTCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGt          <  1:1139570/150‑1 (MQ=255)
                                                                                                                                        aTCAGATAGATGTGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGtt         <  2:1408207/150‑1 (MQ=255)
                                                                                                                                        aTCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGtt         >  1:1517381/1‑150 (MQ=255)
                                                                                                                                           agatagatGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAAt      <  1:667831/150‑1 (MQ=255)
                                                                                                                                               agatGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTggg  <  1:1430980/150‑1 (MQ=255)
                                                                                                                                                    cgcCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCAtc                  >  2:1079899/1‑129 (MQ=255)
                                                                                                                                                    cgcCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCAtc                  <  1:1079899/129‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATCACGCCCGATATCCAGCAAATTGTGATAGCGCACACCAGAAGGACAGGTGGTTTCACAATTACGGCAAGTGAGGCAGCGATCGAGATGCTCCTGTGTTTTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGG  >  NC_000913/3125160‑3125452

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