Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 741,045 A→C 9.6% intergenic (+120/+30) 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,0450AC 9.6% 222.5 / 14.0 84intergenic (+120/+30)phr/dtpDdeoxyribodipyrimidine photolyase, FAD‑binding/dipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base A (37/38);  new base C (3/5);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

CGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAA  >  NC_000913/740897‑741190
                                                                                                                                                    |                                                                                                                                                 
cGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAg                                                                                                                                                  <  2:490403/150‑1 (MQ=255)
    tGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGACCGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATACGACGCGTCAAGCGTc                                                                                                                                              >  1:1595299/1‑150 (MQ=255)
    tGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGAGTCGTTATCGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCATGCGTc                                                                                                                                              <  2:246015/150‑1 (MQ=38)
       ggcggcGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATGCGACGCTTGACGCGTctta                                                                                                                                           >  1:1412658/1‑147 (MQ=25)
         cggcgCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACTTAGGCCTGATAAGACGCTTGACGCGTcttatc                                                                                                                                         <  2:332721/150‑6 (MQ=25)
          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:157162/1‑150 (MQ=255)
                                       aGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTcca                                                                                                           >  1:399543/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:973927/150‑1 (MQ=255)
                                                                                                                            aCGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGt                      <  1:1717280/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)
                                                                                                                                                    |                                                                                                                                                 
CGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAA  >  NC_000913/740897‑741190

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: