Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 741,051 G→T 7.3% intergenic (+126/+24) phr → / ← dtpD deoxyribodipyrimidine photolyase, FAD‑binding/dipeptide and tripeptide permease D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913741,0510GT 7.3% 224.5 / 8.1 82intergenic (+126/+24)phr/dtpDdeoxyribodipyrimidine photolyase, FAD‑binding/dipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base G (38/38);  new base T (2/4);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

GGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAATATATGCATT  >  NC_000913/740907‑741200
                                                                                                                                                |                                                                                                                                                     
ggcgCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCTGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCa                                                                                                                                                  <  2:839652/150‑1 (MQ=255)
  cgcgGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAgg                                                                                                                                                <  2:924491/150‑1 (MQ=255)
   gcgGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATGCGACGCTTGACGCGTCTTATCAGGc                                                                                                                                               <  2:1412658/150‑1 (MQ=25)
   gcgGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATGCGACGCTTGACGCGTCTTATCAGGc                                                                                                                                               <  2:1422930/150‑1 (MQ=25)
    cgGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTTCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                                                                                                              <  2:465858/150‑2 (MQ=25)
    cgGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                                                                                                              <  1:1198188/150‑2 (MQ=25)
        aGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCgg                                                                                                                                          >  1:909070/1‑150 (MQ=255)
          gggAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTg                                                                                                                                        >  1:1487673/1‑150 (MQ=255)
            gAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCt                                                                                                                                      <  1:998964/150‑1 (MQ=255)
             aaaTAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTg                                                                                                                                     <  2:323756/150‑1 (MQ=255)
             aaaTAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTg                                                                                                                                     <  2:718673/150‑1 (MQ=255)
              aaTAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGTCGTGAACGCCTTATCCAGCCTACAAATTTGTTCAAAATCAATAAATTGCAGGAAATACGTAGGCCTGTTAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTgg                                                                                                                                    <  1:1517876/150‑1 (MQ=255)
                tAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACACATTTGTGCAAATTCAATAAATTGCAGGAAATCCGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGATGGca                                                                                                                                  >  1:1307343/1‑149 (MQ=255)
                tAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGcc                                                                                                                                  <  2:145070/150‑1 (MQ=255)
                tAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGcc                                                                                                                                  <  2:129750/150‑1 (MQ=255)
                   cAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTAt                                                                                                                               >  2:1393862/1‑150 (MQ=255)
                       cgcgTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAg                                                                                                                           >  2:1717280/1‑150 (MQ=255)
                       cgcgTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAg                                                                                                                           >  2:825413/1‑150 (MQ=255)
                       cgcgTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAg                                                                                                                           >  2:973927/1‑150 (MQ=255)
                       cgcgTCAGGAAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACTTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAg                                                                                                                           <  2:1210025/150‑1 (MQ=255)
                        gcgTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGa                                                                                                                          <  1:1360405/150‑1 (MQ=255)
                            cAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTcc                                                                                                                      <  1:1230072/150‑1 (MQ=255)
                             aGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTcca                                                                                                                     >  1:399543/1‑150 (MQ=255)
                             aGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTcca                                                                                                                     >  1:157162/1‑150 (MQ=255)
                                cAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATAAGACTccagcgc                                                                                                                 >  1:693613/1‑150 (MQ=255)
                                cAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACCCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcg                                                                                                                  <  1:169850/150‑1 (MQ=255)
                                  aTCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcgcc                                                                                                                >  2:181564/1‑150 (MQ=255)
                                    cGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCCGGCATCGGTGCTGGCCTATTAAGACTccagcgccag                                                                                                              >  1:513813/1‑150 (MQ=255)
                                           aGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCgg                                                                                                       >  1:189245/1‑150 (MQ=255)
                                                     tGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAAtt                                                                                             >  1:949787/1‑150 (MQ=255)
                                                     tGCCGCCTGAACGCCTTCTCCAGCCTACAACTTTGTGCAAATTCAATACATGGCAGGAAATACGTCGGCCTGATCAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGCCTCCAGCGCCAGCGCGCGGTTTCTGAAtt                                                                                             <  2:737337/150‑1 (MQ=255)
                                                       cggcgTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCATTAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTc                                                                                           >  2:649435/1‑150 (MQ=255)
                                                        ggcgTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCa                                                                                          >  2:950405/1‑150 (MQ=255)
                                                        ggcgTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCa                                                                                          <  2:1256716/150‑1 (MQ=255)
                                                         gcgTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCCAGCGTCGCATCAGGCATCGCTGCTGGCATATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAg                                                                                         >  1:1058588/1‑150 (MQ=255)
                                                         gcgTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAg                                                                                         <  1:1637664/150‑1 (MQ=255)
                                                           gTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAGTTTCAgcg                                                                                       >  1:1480316/1‑150 (MQ=255)
                                                            tGAACGCCTTATCCAGCCTAAAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATATGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAgcgc                                                                                      <  1:1087797/150‑1 (MQ=255)
                                                             gAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGACAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGcc                                                                                     >  1:1346767/1‑150 (MQ=255)
                                                                 gCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGAt                                                                                 <  1:1227313/150‑1 (MQ=255)
                                                                     tATCCATCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTATGCCTGATAAGACGCGGCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAg                                                                             <  1:451286/150‑1 (MQ=255)
                                                                     tATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAg                                                                             <  2:1487673/150‑1 (MQ=255)
                                                                        ccAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCa                                                                          >  1:204816/1‑150 (MQ=255)
                                                                               aCAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAatcata                                                                   >  2:629411/1‑149 (MQ=255)
                                                                                cAAATTTGTGAAAATTCAATAAATTTCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAAtcatca                                                                  <  1:1555817/150‑1 (MQ=255)
                                                                                  aaTTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGt                                                                >  2:129982/1‑150 (MQ=255)
                                                                                    tttGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCTGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTAc                                                              >  2:1226988/1‑150 (MQ=255)
                                                                                    tttGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTAc                                                              >  1:945288/1‑150 (MQ=255)
                                                                                      tgtgCAAATTCAATAAATTGCAGGAAATACGTAGCCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGTCCTATTAAGAATCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCa                                                            <  1:1569079/150‑1 (MQ=255)
                                                                                       gtgCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAATTCATCAGTACCAg                                                           >  1:718308/1‑150 (MQ=255)
                                                                                        tgCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGACAAATCATCAGTACCAGt                                                          >  2:1280075/1‑150 (MQ=255)
                                                                                            atttCTATTAATTGAATGATATACGTAGGCCTGATAAGTCGCGTCAAGGGTCGCATCAGGCTTCGGTGCTGGCGTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACGAGTCCTa                                                      <  1:1663087/148‑1 (MQ=255)
                                                                                            aaTTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTa                                                      >  2:669868/1‑150 (MQ=255)
                                                                                               tCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTacac                                                   <  2:399543/150‑1 (MQ=255)
                                                                                               tCAAAAAAATGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTacac                                                   <  2:604426/150‑1 (MQ=255)
                                                                                                cAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTacaca                                                  <  1:1393862/150‑1 (MQ=255)
                                                                                                 aaTAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTacacac                                                 >  2:696478/1‑150 (MQ=255)
                                                                                                  aTAAATTGCAGGACATACGTAGGCCTGATAAGACGCGTCAAGCGGCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGGCCTACACACg                                                <  1:63217/150‑1 (MQ=255)
                                                                                                  aTAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACg                                                <  2:575907/150‑1 (MQ=255)
                                                                                                    aaaTTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGcc                                              >  2:1376719/1‑150 (MQ=255)
                                                                                                                atttgTAGGCTTGATAAGGCGCGTCAAGCTGCGCATCTGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAg                                  <  2:726714/146‑1 (MQ=255)
                                                                                                                  aCGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGt                                <  1:1717280/150‑1 (MQ=255)
                                                                                                                  aCGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGt                                <  1:973927/150‑1 (MQ=255)
                                                                                                                   cGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTa                               >  2:370050/1‑150 (MQ=255)
                                                                                                                     tAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGACAGTGCGCCCCAGGTAAt                             >  1:823164/1‑150 (MQ=255)
                                                                                                                       tgCCTGATGCGACGCTTGACGCGTCTTATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAttt                           >  1:580107/2‑150 (MQ=255)
                                                                                                                       tgCCTGATGCGACGCTTGACGCGTCTTATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAttt                           >  1:587470/2‑150 (MQ=255)
                                                                                                                          cTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGAt                        >  1:540662/1‑150 (MQ=255)
                                                                                                                          cTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGAt                        >  1:642604/1‑150 (MQ=255)
                                                                                                                           tGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATc                       <  2:112404/150‑1 (MQ=255)
                                                                                                                           tGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATc                       <  1:818548/150‑1 (MQ=255)
                                                                                                                           tGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATAAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCACGTAATTTGATc                       <  2:1395681/150‑1 (MQ=255)
                                                                                                                            taTAAGACTCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCa                      <  1:881946/149‑1 (MQ=255)
                                                                                                                             aTAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCaa                     <  2:899107/150‑1 (MQ=255)
                                                                                                                             aTAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCaa                     >  1:1078691/1‑150 (MQ=255)
                                                                                                                                aGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATa                  >  2:1116844/1‑150 (MQ=255)
                                                                                                                                     cgTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTCCCAGTCCTACCCACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCa             >  1:908856/1‑150 (MQ=255)
                                                                                                                                      gTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCaa            <  1:601518/150‑1 (MQ=255)
                                                                                                                                      gTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCaa            <  2:946902/150‑1 (MQ=255)
                                                                                                                                            cGTCGCATCAGGCAGCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATAATAAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAATATATg      <  2:1710329/150‑1 (MQ=255)
                                                                                                                                             gTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAATATATGc     >  2:283164/1‑150 (MQ=255)
                                                                                                                                             gTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAATATATGc     >  2:1530119/1‑150 (MQ=255)
                                                                                                                                                gCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAATATATGCAtt  <  1:1326924/150‑1 (MQ=255)
                                                                                                                                                gCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAATATATGCAtt  <  2:1595299/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAATATATGCATT  >  NC_000913/740907‑741200

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