Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,074,814:1 +A 5.3% coding (36/1533 nt) yeeR → CP4‑44 prophage; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,074,8141.A 5.3% 572.0 / 26.0 133coding (36/1533 nt)yeeRCP4‑44 prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base . (73/53);  new base A (5/2);  total (78/55)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

Only 100 of 136 total aligned reads displayed.

CGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCATATTCGC  >  NC_000913/2074668‑2074961
                                                                                                                                                   |                                                                                                                                                   
cGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑g                                                                                                                                                    >  2:616342/1‑148 (MQ=255)
cGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑g                                                                                                                                                    >  2:853604/1‑148 (MQ=255)
   ctctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACTGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCatg                                                                                                                                                >  3:874139/1‑147 (MQ=255)
   ctctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAg                                                                                                                                                 >  1:60735/1‑148 (MQ=255)
     ctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGa                                                                                                                                               >  1:216604/1‑148 (MQ=255)
      tctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGAtt                                                                                                                                             >  4:845000/1‑149 (MQ=255)
             tCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGa                                                                                                                                               <  4:400781/140‑1 (MQ=255)
             tCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAtt                                                                                                                                      >  4:714443/1‑149 (MQ=255)
               ccGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGGTAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttct                                                                                                                                    >  4:561712/1‑149 (MQ=255)
               ccGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttct                                                                                                                                    >  3:446243/1‑149 (MQ=255)
                cGGTCATGATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttct                                                                                                                                    >  2:708013/1‑148 (MQ=255)
                cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttct                                                                                                                                    >  1:913029/1‑148 (MQ=255)
                   tcatgaTGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCg                                                                                                                                 <  1:708013/148‑1 (MQ=255)
                    catcatGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGc                                                                                                                                >  1:311199/1‑148 (MQ=255)
                      tcatGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCtt                                                                                                                             <  4:660309/149‑1 (MQ=255)
                      tcatGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCt                                                                                                                              <  2:673855/148‑1 (MQ=255)
                        atGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCttt                                                                                                                            <  1:410293/148‑1 (MQ=255)
                             cGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCa                                                                                                                      <  3:530181/149‑1 (MQ=255)
                              gggACCCGGACCGCCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCa                                                                                                                      >  2:742506/1‑148 (MQ=255)
                                   ccGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCAt                                                                                                                >  4:32513/1‑149 (MQ=255)
                                   ccGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCAt                                                                                                                >  3:250504/1‑149 (MQ=255)
                                       aCCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAAt                                                                                                            <  4:250504/149‑1 (MQ=255)
                                         cGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATc                                                                                                           >  1:822552/1‑148 (MQ=255)
                                              cAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTg                                                                                                                                           <  1:91251/111‑1 (MQ=255)
                                              cAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTg                                                                                                                                           >  2:91251/1‑111 (MQ=255)
                                                aCGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACgg                                                                                                    <  2:216604/148‑1 (MQ=255)
                                                  ggATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCtt                                                                                                 >  3:100385/1‑149 (MQ=255)
                                                   gATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttcttc                                                                                                                                  <  3:399853/115‑1 (MQ=255)
                                                   gATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTc                                                                                                <  3:309963/149‑1 (MQ=255)
                                                    aTCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTc                                                                                                <  2:311199/148‑1 (MQ=255)
                                                    aTCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCt                                                                                               >  3:197277/1‑149 (MQ=255)
                                                       ttCAACGCCACATTCGCTGGCATTAACAGTAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTgcg                                                                                            >  3:777248/1‑149 (MQ=255)
                                                         cAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTgcgc                                                                                           >  2:22859/1‑148 (MQ=255)
                                                           aCGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTg                                                                                         >  1:537025/1‑148 (MQ=255)
                                                            cGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAAt                                                                                                            <  2:884884/128‑1 (MQ=255)
                                                               cacaTTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAgg                                                                                    >  3:253532/1‑149 (MQ=255)
                                                                 caTTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAggg                                                                                   <  1:616342/148‑1 (MQ=255)
                                                                 caTTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGc                                                                                  <  3:714443/149‑1 (MQ=255)
                                                                  aTTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCt                                                                                 <  4:874139/149‑1 (MQ=255)
                                                                   ttCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCt                                                                                 >  2:13261/1‑148 (MQ=255)
                                                                        tGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATAtt                                                                           >  3:756118/1‑149 (MQ=255)
                                                                          gCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                                                         <  3:845000/149‑1 (MQ=255)
                                                                           cATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                                                         <  2:822552/148‑1 (MQ=255)
                                                                           cATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGACATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCt                                                                        <  3:561712/149‑1 (MQ=255)
                                                                           cATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTACTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCt                                                                        <  4:149415/149‑1 (MQ=255)
                                                                            aTTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCt                                                                                 >  3:266849/1‑139 (MQ=255)
                                                                             ttAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtgt                                                                      >  3:589085/1‑149 (MQ=255)
                                                                             ttAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtgt                                                                      >  4:274492/1‑149 (MQ=255)
                                                                             ttAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtgt                                                                      >  3:369419/1‑149 (MQ=255)
                                                                             ttAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtg                                                                       >  1:70580/1‑148 (MQ=255)
                                                                             ttAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGC‑GATCAGCAGGAt                                                                                                                                              >  2:54124/1‑77 (MQ=255)
                                                                             ttAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGC‑GATCAGCAGGAt                                                                                                                                              <  1:54124/77‑1 (MQ=255)
                                                                              tAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtgt                                                                      >  1:323003/1‑148 (MQ=255)
                                                                                aCAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTc                                                                   >  4:85307/1‑149 (MQ=255)
                                                                                aCAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTc                                                                   <  4:777248/149‑1 (MQ=255)
                                                                                   aTAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                                                         <  1:433573/140‑1 (MQ=255)
                                                                                   aTAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                                                         >  2:433573/1‑140 (MQ=255)
                                                                                       cATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCa                                                                                                                      <  3:430762/91‑1 (MQ=255)
                                                                                       cATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCa                                                                                                                      >  4:430762/1‑91 (MQ=255)
                                                                                        aTGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTc                                                                                                >  1:682785/1‑112 (MQ=255)
                                                                                         tGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggc                                                           <  2:318720/148‑1 (MQ=255)
                                                                                         tGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggc                                                           >  2:837590/1‑148 (MQ=255)
                                                                                          gATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcc                                                          >  1:804602/1‑148 (MQ=255)
                                                                                            tatTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtgt                                                                      >  3:472132/1‑134 (MQ=255)
                                                                                            tatTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtgt                                                                      <  4:472132/134‑1 (MQ=255)
                                                                                                  tcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTAc                                                 >  4:426943/1‑149 (MQ=255)
                                                                                                        aGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActggctg                                           <  3:638093/149‑1 (MQ=255)
                                                                                                        aGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActggct                                            <  1:853604/148‑1 (MQ=255)
                                                                                                          tGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTTCTGTGTCTGTTcggcccggc                                                      <  4:756402/136‑1 (MQ=255)
                                                                                                          tGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcccggc                                                      >  3:756402/1‑136 (MQ=255)
                                                                                                           gACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTggctggct                                        >  4:657054/1‑149 (MQ=255)
                                                                                                            aCTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCta                                       >  4:775952/1‑149 (MQ=255)
                                                                                                              tATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcc                                                          <  3:807754/128‑1 (MQ=255)
                                                                                                              tATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcc                                                          >  4:807754/1‑128 (MQ=255)
                                                                                                               aTGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTc                                                                                                <  2:30594/89‑1 (MQ=255)
                                                                                                               aTGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTc                                                                                                >  1:30594/1‑89 (MQ=255)
                                                                                                               aTGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCtata                                     >  1:810412/1‑148 (MQ=255)
                                                                                                                        aTAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAAc                            >  1:691182/1‑148 (MQ=255)
                                                                                                                        aTAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACg                           >  4:866958/1‑149 (MQ=255)
                                                                                                                          aGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGc                          >  1:674674/1‑148 (MQ=255)
                                                                                                                             cGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCACCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATa                       >  2:323127/1‑148 (MQ=255)
                                                                                                                               gcgcgcTGATCCTGCTGATC‑GCAGGATTTGCAATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAAc                     <  1:429630/148‑1 (MQ=255)
                                                                                                                                    cTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGgt                >  1:863352/1‑148 (MQ=255)
                                                                                                                                     aaaTCCTGC‑GATCAGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGgttg              >  3:416535/3‑149 (MQ=255)
                                                                                                                                     aaaTCCTGC‑GATCAGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGgttg              >  3:493703/3‑149 (MQ=255)
                                                                                                                                      aaTCCTGC‑GATCAGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGgttgt             >  4:369059/2‑149 (MQ=255)
                                                                                                                                       aTCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGgttgt             <  2:60735/148‑1 (MQ=255)
                                                                                                                                       aTCCTGC‑GATCAGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGgt                <  2:696686/145‑1 (MQ=255)
                                                                                                                                         ccTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTcc          <  4:253532/149‑1 (MQ=255)
                                                                                                                                            gctgATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                                                         <  2:549376/83‑1 (MQ=255)
                                                                                                                                            gctgATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                                                         >  1:549376/1‑83 (MQ=255)
                                                                                                                                            gctgATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCata       >  3:905107/1‑149 (MQ=255)
                                                                                                                                            gctgATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCata       <  3:866958/149‑1 (MQ=255)
                                                                                                                                            gctgATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCat        >  1:199556/1‑148 (MQ=255)
                                                                                                                                              tgATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCATAtt     >  3:194011/1‑149 (MQ=255)
                                                                                                                                               gATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCATAtt     >  2:333602/1‑148 (MQ=255)
                                                                                                                                                aTC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCATATTCg   >  3:626370/1‑149 (MQ=255)
                                                                                                                                                 tC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCta                                       <  2:15534/112‑1 (MQ=255)
                                                                                                                                                 tC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCta                                       >  1:15534/1‑112 (MQ=255)
                                                                                                                                                 tC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCATATTCg   >  1:822322/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCATATTCGC  >  NC_000913/2074668‑2074961

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