Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 741,043 C→G 9.2% intergenic (+118/+32) 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,0430CG 9.2% 271.0 / 16.3 87intergenic (+118/+32)phr/dtpDdeoxyribodipyrimidine photolyase, FAD‑binding/dipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base C (37/42);  new base G (3/5);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

CAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAA  >  NC_000913/740894‑741190
                                                                                                                                                     |                                                                                                                                                   
cAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCTGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTc                                                                                                                                                     >  1:839652/1‑150 (MQ=255)
cAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTc                                                                                                                                                     >  2:67433/1‑150 (MQ=255)
 aGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCa                                                                                                                                                    <  2:1166711/150‑1 (MQ=255)
 aGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCa                                                                                                                                                    <  2:1314231/150‑1 (MQ=255)
 aGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCa                                                                                                                                                    <  2:1670728/150‑1 (MQ=255)
   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:1422930/150‑1 (MQ=25)
                gcgGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATGCGACGCTTGACGCGTCTTATCAGGc                                                                                                                                     <  2:1412658/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:129750/150‑1 (MQ=255)
                             tAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGcc                                                                                                                        <  2:145070/150‑1 (MQ=255)
                                cAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTAt                                                                                                                     >  2:1393862/1‑150 (MQ=255)
                                    cgcgTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAg                                                                                                                 >  2:825413/1‑150 (MQ=255)
                                    cgcgTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAg                                                                                                                 >  2:1717280/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: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:587470/2‑150 (MQ=255)
                                                                                                                                    tgCCTGATGCGACGCTTGACGCGTCTTATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAAttt                 >  1:580107/2‑150 (MQ=255)
                                                                                                                                       cTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGAt              >  1:642604/1‑150 (MQ=255)
                                                                                                                                       cTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGAt              >  1:540662/1‑150 (MQ=255)
                                                                                                                                        tGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATc             <  1:818548/150‑1 (MQ=255)
                                                                                                                                        tGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATc             <  2:112404/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  <  2:946902/150‑1 (MQ=255)
                                                                                                                                                   gTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCaa  <  1:601518/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCATCAGTACCAGTCCTACACACGCCAGTGCGCCCCAGGTAATTTGATCAAATACTTCAA  >  NC_000913/740894‑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: