Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,889 T→A 10.9% intergenic (+33/‑50) ygeI → / → pbl protein YgeI/transglycosylase domain‑containing protein Pbl

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,993,8890TA10.9% 586.5 / 43.1 200intergenic (+33/‑50)ygeI/pblprotein YgeI/transglycosylase domain‑containing protein Pbl
Reads supporting (aligned to +/- strand):  ref base T (114/64);  new base A (13/9);  total (127/73)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

GCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGA  >  NC_000913/2993740‑2994037
                                                                                                                                                     |                                                                                                                                                    
     tGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                                                                                                                                                 <  1:241944/150‑1 (MQ=255)
      gATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataa                                                                                                                                                >  2:1232490/1‑150 (MQ=255)
          cAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGt                                                                                                                                            <  2:21129/150‑1 (MQ=255)
           aGATGAAATACAAGCATTTATAGTTTCAACATGTATTGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                                                                           <  1:740021/150‑1 (MQ=255)
                 aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                                                                     <  2:1355365/150‑1 (MQ=255)
                  aTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGtt                                                                                                                                    <  1:1530141/150‑1 (MQ=255)
                   tACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTc                                                                                                                                   >  2:177539/1‑150 (MQ=255)
                    aCAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc                                                                                                                                  <  1:427474/150‑1 (MQ=255)
                    aCAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc                                                                                                                                  <  2:1375207/150‑1 (MQ=255)
                    aCAACCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                                                                                  <  1:839093/150‑3 (MQ=255)
                     cAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatg                                                                                                                                 <  2:832701/150‑4 (MQ=255)
                        gCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                                                                              >  1:500125/1‑150 (MQ=255)
                         cATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttca                                                                                                                             <  2:169105/150‑1 (MQ=255)
                            ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattc                                                                                                                          <  1:527479/150‑1 (MQ=255)
                             tataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattca                                                                                                                         >  1:301846/1‑150 (MQ=255)
                                aGTTTCAACTTGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                                                                      >  1:1478349/1‑150 (MQ=255)
                                 gTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                                                                     >  2:1286692/1‑150 (MQ=255)
                                 gTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                                                                     <  1:633012/150‑1 (MQ=255)
                                     cAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                                                                                 <  1:718012/150‑1 (MQ=255)
                                       aCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacac                                                                                                               >  2:1187119/1‑150 (MQ=255)
                                        cATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacaca                                                                                                              <  2:563704/150‑1 (MQ=255)
                                          tggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATg                                                                                                            >  1:694629/1‑150 (MQ=255)
                                          tggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATg                                                                                                            <  2:849619/150‑1 (MQ=255)
                                               ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaaa                                                                                                       >  2:305778/1‑150 (MQ=255)
                                                  gcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                                                                    >  1:352736/1‑150 (MQ=255)
                                                   ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAatat                                                                                                   >  2:622676/1‑150 (MQ=255)
                                                       ttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTACTTTCATTAAATCCTCACACATGAAAAAATATgcta                                                                                               >  2:250345/1‑148 (MQ=255)
                                                        tCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCt                                                                                              >  2:154680/1‑150 (MQ=255)
                                                              tGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaa                                                                                        >  1:673281/1‑150 (MQ=255)
                                                                 tatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGATGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                                                                     >  1:685688/1‑150 (MQ=255)
                                                                 tatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                                                                     >  1:323701/1‑150 (MQ=255)
                                                                   tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGtt                                                                                   >  2:75916/1‑150 (MQ=255)
                                                                   tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGtt                                                                                   >  2:1220885/1‑150 (MQ=255)
                                                                        tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattc                                                                                                                          <  2:648260/106‑1 (MQ=255)
                                                                           gaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCtt                                                                           >  1:105373/1‑150 (MQ=255)
                                                                           gaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCtt                                                                           <  1:154680/150‑1 (MQ=255)
                                                                             agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTaa                                                                         >  2:1497144/1‑150 (MQ=255)
                                                                             agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTaa                                                                         >  2:459578/1‑150 (MQ=255)
                                                                             agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTaa                                                                         >  2:886452/1‑150 (MQ=255)
                                                                              gatgatCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                                                                        <  2:391066/150‑1 (MQ=255)
                                                                              gaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                                                                        <  1:134309/150‑1 (MQ=255)
                                                                              gaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                                                                        <  1:797531/150‑1 (MQ=255)
                                                                              gaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                                                                        >  2:797531/1‑150 (MQ=255)
                                                                                  aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                                                                    >  1:35269/1‑150 (MQ=255)
                                                                                  aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                                                                    >  1:767144/1‑150 (MQ=255)
                                                                                    catATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGAGCTCACCAAAAGTTAAATTATTAACCAACta                                                                  >  2:663098/3‑149 (MQ=255)
                                                                                        tctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg                                                              <  1:1241037/150‑1 (MQ=255)
                                                                                         ctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt                                                             >  1:746689/1‑150 (MQ=255)
                                                                                           cACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCGCGttt                                                           >  2:1322310/1‑150 (MQ=255)
                                                                                             ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg                                                         >  2:950407/1‑150 (MQ=255)
                                                                                              attACTATAAAATTTAGTTGTAGACGATTTACAAAATAAATCATAACTAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGt                                                        <  1:1356430/149‑1 (MQ=255)
                                                                                                tACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtat                                                      >  2:26746/1‑150 (MQ=255)
                                                                                                tACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtat                                                      <  1:622676/150‑1 (MQ=255)
                                                                                                 aCTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata                                                     <  1:912502/150‑1 (MQ=255)
                                                                                                  cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTGCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAAGTTTTGTatat                                                    >  2:296093/1‑150 (MQ=255)
                                                                                                  cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTatat                                                    >  2:1172861/1‑150 (MQ=255)
                                                                                                   tataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatt                                                   >  2:1413273/1‑150 (MQ=255)
                                                                                                   aataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatt                                                   <  2:1416554/149‑1 (MQ=255)
                                                                                                    ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtattt                                                  >  2:1121687/1‑150 (MQ=255)
                                                                                                    ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtattt                                                  >  1:817414/1‑150 (MQ=255)
                                                                                                    ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtattt                                                  >  1:679271/1‑150 (MQ=255)
                                                                                                      aaaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatttat                                                >  1:1277648/1‑150 (MQ=255)
                                                                                                       aaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatttatt                                               >  2:1304227/1‑150 (MQ=255)
                                                                                                            gattGGTAGACGATTGAAAAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACaaa                                          >  2:1049208/4‑150 (MQ=255)
                                                                                                             aGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACaaaa                                         <  2:1061329/150‑1 (MQ=255)
                                                                                                              gTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACaaaaa                                        >  2:397440/1‑150 (MQ=255)
                                                                                                              gTGGGAGACGATTGACCAAATAAATCATAACGAGATACTTTAATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTGTTGTATATTTATTACaaaaa                                        <  2:806750/150‑1 (MQ=255)
                                                                                                                ggTAGACGATTTACCAAATAAATCATAACGTGATACTTTTATAATTGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATg                                      <  1:584221/150‑1 (MQ=255)
                                                                                                                ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATg                                      >  1:1388759/1‑150 (MQ=255)
                                                                                                                 gTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGt                                     >  1:216288/1‑150 (MQ=255)
                                                                                                                 gTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGt                                     >  2:615264/1‑150 (MQ=255)
                                                                                                                   aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCt                                   >  2:586556/1‑150 (MQ=255)
                                                                                                                   aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCt                                   >  1:775175/1‑150 (MQ=255)
                                                                                                                      cGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTctc                                                                                           <  2:509187/91‑3 (MQ=255)
                                                                                                                            aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGc                          >  1:124290/1‑150 (MQ=255)
                                                                                                                            aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGc                          >  2:356038/1‑150 (MQ=255)
                                                                                                                            aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGc                          <  2:727419/150‑1 (MQ=255)
                                                                                                                               aaataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGAGCTAACCAAAAGTTAAATTATTAACAATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCaaa                       >  1:1450713/1‑149 (MQ=255)
                                                                                                                                    gaaaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGAGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTAAAAAAATGTATTTCGCCAGCAACTCTTg                  >  1:641442/7‑150 (MQ=39)
                                                                                                                                     atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATGCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGc                 <  2:1265722/150‑1 (MQ=255)
                                                                                                                                     atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGc                 >  1:170020/1‑150 (MQ=255)
                                                                                                                                      tCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGct                <  2:1478349/150‑1 (MQ=255)
                                                                                                                                       cATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGctc               <  2:760604/150‑1 (MQ=255)
                                                                                                                                        aTAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGctct              <  1:795814/150‑1 (MQ=255)
                                                                                                                                        aTAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGctct              <  1:939859/150‑1 (MQ=255)
                                                                                                                                         gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCTCCAGCAACTCTTGCTCTa             >  1:1034135/2‑150 (MQ=255)
                                                                                                                                         gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTa             >  1:1141068/2‑150 (MQ=255)
                                                                                                                                         gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTa             >  1:1332993/2‑150 (MQ=255)
                                                                                                                                         caaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCATCTCTTGCTCTa             >  2:1441948/2‑150 (MQ=255)
                                                                                                                                           aCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAAt           >  2:1533458/1‑150 (MQ=255)
                                                                                                                                           aCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAAt           >  2:795166/1‑150 (MQ=255)
                                                                                                                                           aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAAt           >  1:131676/1‑150 (MQ=255)
                                                                                                                                           aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAAt           >  1:1366423/1‑150 (MQ=255)
                                                                                                                                            cGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATg          <  1:184490/150‑1 (MQ=255)
                                                                                                                                             gagaTACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGa         >  1:142011/1‑150 (MQ=255)
                                                                                                                                              agaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGaa        <  1:1533458/150‑1 (MQ=255)
                                                                                                                                                 tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCa     >  2:424367/1‑150 (MQ=255)
                                                                                                                                                  aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAg    <  1:75916/150‑1 (MQ=255)
                                                                                                                                                   cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAgg   >  1:1069027/1‑150 (MQ=255)
                                                                                                                                                   cTTTTATAATAGTATCTCGTTCATTACATTCAAACCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAACATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCCCCAGCAACTCTTGCTCTAAAGAAGCAgg   >  2:691279/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGA  >  NC_000913/2993740‑2994037

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