Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,894,852 T→G 6.1% intergenic (+81/+67) ygcP → / ← ygcQ putative anti‑terminator regulatory protein/putative flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,894,8520TG 6.1% 413.5 / 16.4 131intergenic (+81/+67)ygcP/ygcQputative anti‑terminator regulatory protein/putative flavoprotein
Reads supporting (aligned to +/- strand):  ref base T (60/63);  new base G (2/6);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 133 total aligned reads displayed.

AACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTG  >  NC_000913/2894703‑2895001
                                                                                                                                                     |                                                                                                                                                     
aaCGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCt                                                                                                                                                       <  2:1260824/150‑1 (MQ=255)
    cgGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAc                                                                                                                                                   <  1:510210/150‑1 (MQ=255)
     gggTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcg                                                                                                                                                  <  2:1488553/150‑1 (MQ=255)
      ggTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                                                                                                                 >  1:252048/1‑150 (MQ=255)
         gtCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTc                                                                                                                                              <  1:1074314/150‑1 (MQ=255)
           cGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCtt                                                                                                                                            <  1:1106379/150‑1 (MQ=255)
           cGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCtt                                                                                                                                            >  2:1278872/1‑150 (MQ=255)
            gTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTa                                                                                                                                           >  1:1381767/1‑150 (MQ=255)
              ggCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTTGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                                                                         <  2:1495206/150‑1 (MQ=255)
                cgcTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAg                                                                                                                                       <  1:990521/150‑1 (MQ=255)
                 gcattCCACCACGAATAACCGCGTCTGGACGTTATCGAAAAAAATACTTTGACGGGATAATAGTAAACAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                                                                      >  2:1477472/4‑150 (MQ=255)
                 gcTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                                                                      <  1:1358116/150‑1 (MQ=255)
                  cTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGc                                                                                                                                     <  1:1492286/150‑1 (MQ=255)
                  cTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGc                                                                                                                                     <  1:1337592/150‑1 (MQ=255)
                    ttCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgc                                                                                                                                   <  2:366381/150‑3 (MQ=255)
                    ttCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgc                                                                                                                                   <  1:1256286/150‑3 (MQ=255)
                    ttCCAACACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgc                                                                                                                                   <  1:429583/150‑3 (MQ=255)
                      ccaccaCGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                                                                                                                 >  1:1563425/1‑150 (MQ=255)
                       caccacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGTGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACa                                                                                                                                >  1:1003230/1‑150 (MQ=255)
                       caccacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACa                                                                                                                                >  2:796695/1‑150 (MQ=255)
                       caccacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACa                                                                                                                                >  2:601279/1‑150 (MQ=255)
                        accacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACaa                                                                                                                               >  1:232172/1‑150 (MQ=255)
                         ccacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAg                                                                                                                              >  2:514663/1‑150 (MQ=255)
                         ccacGAATACCGGGGTCTGGACGTTAGCGAAAAAATAACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAg                                                                                                                              <  1:1069236/150‑1 (MQ=255)
                          cacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGTTAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGt                                                                                                                             >  2:682016/1‑150 (MQ=255)
                           acGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTc                                                                                                                            >  1:1290480/1‑150 (MQ=255)
                                tACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAg                                                                                                                       <  2:1254045/150‑1 (MQ=255)
                                  ccGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGc                                                                                                                     <  2:1097607/150‑1 (MQ=255)
                                   cGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGTCTACAAGTCCCCAGGcc                                                                                                                    <  1:620572/150‑1 (MQ=255)
                                   cGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGcc                                                                                                                    >  2:1443786/1‑150 (MQ=255)
                                   cGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGcc                                                                                                                    >  2:1449700/1‑150 (MQ=255)
                                         cTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAAcc                                                                                                              <  1:1566309/150‑1 (MQ=255)
                                           ggACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGt                                                                                                            >  1:1239614/1‑150 (MQ=255)
                                            gACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTa                                                                                                           <  2:274748/150‑1 (MQ=255)
                                              cGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAgg                                                                                                         >  1:878879/1‑150 (MQ=255)
                                               gTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGc                                                                                                        >  2:1631814/1‑150 (MQ=255)
                                               gTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGc                                                                                                        <  2:299642/150‑1 (MQ=255)
                                                ttAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCa                                                                                                       >  1:977817/1‑150 (MQ=255)
                                                 tAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAg                                                                                                      >  1:456977/1‑150 (MQ=255)
                                                   gCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGa                                                                                                    >  1:117243/1‑150 (MQ=255)
                                                   gCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGa                                                                                                    >  1:736882/1‑150 (MQ=255)
                                                    cGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGat                                                                                                   >  2:1312116/1‑150 (MQ=255)
                                                        aaaaTTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGatcatc                                                                                               >  2:602249/1‑150 (MQ=255)
                                                            ttACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGc                                                                                           <  2:1088915/150‑1 (MQ=255)
                                                                  ttaCGGGATAATCGTTACCTATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACg                                                                                     <  2:602868/148‑1 (MQ=255)
                                                                       ggATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGTCATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTAtt                                                                                <  2:1437605/150‑1 (MQ=255)
                                                                       ggATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                                                                      >  1:1088857/1‑96 (MQ=255)
                                                                       ggATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                                                                      <  2:1088857/96‑1 (MQ=255)
                                                                       ggATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTAtt                                                                                <  1:1287529/150‑1 (MQ=255)
                                                                        gATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCAGATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCGCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTg                                                                               <  2:399902/150‑1 (MQ=255)
                                                                          tAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGc                                                                             <  2:1378843/150‑1 (MQ=255)
                                                                           aaTCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCa                                                                            >  2:928199/1‑150 (MQ=255)
                                                                            aTCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCaa                                                                           >  1:1412061/1‑150 (MQ=255)
                                                                            aTCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGATATTGGCaa                                                                           >  2:344505/1‑150 (MQ=255)
                                                                                    cAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTg                                                                   <  1:1292257/150‑1 (MQ=255)
                                                                                    cAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTg                                                                   <  2:126884/150‑1 (MQ=255)
                                                                                    cAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTg                                                                   <  2:736882/150‑1 (MQ=255)
                                                                                       ttGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCAGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAAt                                                                >  1:138500/1‑150 (MQ=255)
                                                                                        tGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATg                                                               <  1:1291921/150‑1 (MQ=255)
                                                                                           aTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGttt                                                            <  2:1601257/150‑1 (MQ=255)
                                                                                            ttttGTTTGTTTTGTATGGCGGACCGCAAAATCGTCATTCACCGCCTGAATCGTAGCTTGCAGGGTCTTATCATGCCTAAAATGCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTg                                                           <  2:1290480/150‑1 (MQ=255)
                                                                                            tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTg                                                           <  2:696549/150‑1 (MQ=255)
                                                                                            tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTg                                                           <  2:117243/150‑1 (MQ=255)
                                                                                            tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTg                                                           <  2:364073/150‑1 (MQ=255)
                                                                                             ttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGt                                                          <  1:1631814/150‑1 (MQ=255)
                                                                                              tGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTc                                                         <  2:1304839/150‑1 (MQ=255)
                                                                                               ggTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCa                                                        <  1:126327/150‑1 (MQ=255)
                                                                                                 ttggATTTGTAGGCCTCACGCCAAATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTGTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCAcc                                                      <  1:710398/147‑1 (MQ=255)
                                                                                                    gATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGt                                                   >  2:1107171/1‑150 (MQ=255)
                                                                                                      tttGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGc                                                 >  1:136616/1‑150 (MQ=255)
                                                                                                          tAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTcc                                             >  2:796202/1‑150 (MQ=255)
                                                                                                            ggCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAg                                           >  2:843876/1‑150 (MQ=255)
                                                                                                              ccGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGca                                         <  1:1170362/150‑1 (MQ=255)
                                                                                                              ccGCACGCCACATACGAAATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCATGCCTACAAGTCCCCAGTAGAGAACCGTAGGCAGGATCATCCTGCTACACGTTATTGGCAATCTGCGTTAATGTTTGTGACCAGTGGGTCCTGca                                         <  2:465523/150‑1 (MQ=37)
                                                                                                               cGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGcac                                        <  1:682016/150‑1 (MQ=255)
                                                                                                                gCACGCCACATCCGACCTTCAGCGACTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCAGCCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATAGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTccaccacca                                      <  2:1141378/150‑1 (MQ=38)
                                                                                                                    gCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACtt                                   >  1:1333283/1‑150 (MQ=255)
                                                                                                                     ccACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACttt                                  >  2:472720/1‑150 (MQ=255)
                                                                                                                        caTCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCa                               <  1:1337419/150‑1 (MQ=255)
                                                                                                                         aTCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAg                              <  1:602249/150‑1 (MQ=255)
                                                                                                                           ccGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtc                            >  1:1375675/1‑150 (MQ=255)
                                                                                                                            cGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGTAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtca                           >  2:468279/1‑150 (MQ=255)
                                                                                                                            cGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGGGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtca                           <  1:1218653/150‑1 (MQ=255)
                                                                                                                              aCATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtcatc                         <  1:1302737/150‑1 (MQ=255)
                                                                                                                               cATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtcatca                        <  1:1481416/150‑1 (MQ=255)
                                                                                                                                   cAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCa                    >  2:1432254/1‑150 (MQ=255)
                                                                                                                                    aGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCAc                   <  2:463646/150‑1 (MQ=255)
                                                                                                                                     gcgcCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACg                  >  2:955638/1‑150 (MQ=255)
                                                                                                                                     gcgcCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCACCAGCACCACTTTCCAGACATCAACCAcc                  >  2:1070104/1‑149 (MQ=255)
                                                                                                                                     aggcCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACg                  >  2:540876/3‑150 (MQ=255)
                                                                                                                                      ggcCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGc                 >  2:7661/2‑150 (MQ=255)
                                                                                                                                          tGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAAc             <  2:977817/150‑1 (MQ=255)
                                                                                                                                           gATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAAt                                                                          <  1:1461239/88‑1 (MQ=255)
                                                                                                                                           gATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAAt                                                                          >  2:1461239/1‑88 (MQ=255)
                                                                                                                                           gATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACa            <  2:878879/150‑1 (MQ=255)
                                                                                                                                              gCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAAc             <  1:1545339/146‑1 (MQ=255)
                                                                                                                                              gCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAAc             >  2:1545339/1‑146 (MQ=255)
                                                                                                                                               cGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTg        >  2:185749/1‑150 (MQ=255)
                                                                                                                                                    cTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCtgt   >  2:945296/1‑150 (MQ=255)
                                                                                                                                                     ttGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCtgtg  >  1:762363/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTG  >  NC_000913/2894703‑2895001

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