Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,792,767 A→T 8.2% intergenic (‑216/+59) kbl ← / ← yibB glycine C‑acetyltransferase/YibB family protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,792,7670AT 8.2% 604.9 / 31.6 171intergenic (‑216/+59)kbl/yibBglycine C‑acetyltransferase/YibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (78/79);  new base T (7/7);  total (85/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 186 total aligned reads displayed.

GCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTT  >  NC_000913/3792619‑3792913
                                                                                                                                                    |                                                                                                                                                  
         cTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtt                                                                                                                                           <  2:1249529/149‑1 (MQ=255)
          ttCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACAACTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                          <  1:1455494/149‑1 (MQ=255)
               cAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCt                                                                                                                                     >  4:837571/1‑149 (MQ=255)
                     tCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTCTTCCTATCTCa                                                                                                                               <  3:837571/149‑1 (MQ=255)
                         gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  4:295250/125‑1 (MQ=255)
                         gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  4:1437046/125‑1 (MQ=255)
                         gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  4:1070572/125‑1 (MQ=255)
                         gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  3:1437046/1‑125 (MQ=255)
                         gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  3:1070572/1‑125 (MQ=255)
                         gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTg                                                                                                                           >  1:1025933/1‑149 (MQ=255)
                               taaCGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTc                                                                                                                                       >  2:638921/1‑131 (MQ=255)
                                    tgtAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  4:714861/1‑114 (MQ=255)
                                    tgtAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  3:714861/114‑1 (MQ=255)
                                                  tctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATgg                                                                                                  >  1:1077247/1‑149 (MQ=255)
                                                    tcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                          >  2:895130/1‑107 (MQ=255)
                                                       acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGa                                                                                             <  1:1446620/149‑1 (MQ=255)
                                                         acCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcat                                                                                                                                               >  2:767206/1‑97 (MQ=255)
                                                         acCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcat                                                                                                                                               <  1:767206/97‑1 (MQ=255)
                                                          ccTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAAc                                                                                          <  3:318209/149‑1 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                           >  1:883281/1‑99 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatca                                                                                                                                             <  2:18572/97‑1 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatca                                                                                                                                             >  3:271174/1‑97 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAttttt                                                                                                                                        >  2:710969/1‑102 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATc                                                                                                                      >  2:624290/1‑120 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATc                                                                                                                      <  1:624290/120‑1 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAAcc                                                                                         >  3:517097/1‑149 (MQ=255)
                                                            tACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                                >  4:654488/1‑93 (MQ=255)
                                                            tACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCt                                                                                        <  3:886550/149‑1 (MQ=255)
                                                               aGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                                <  3:743230/90‑1 (MQ=255)
                                                               aGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                                >  4:743230/1‑90 (MQ=255)
                                                                gATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaaa                                                                                           <  1:406176/142‑1 (MQ=255)
                                                                gATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaaa                                                                                           >  2:406176/1‑142 (MQ=255)
                                                                     tagAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGt                                                                                                                         >  1:567825/1‑107 (MQ=255)
                                                                       gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAttttt                                                                                                                                        >  3:26749/1‑90 (MQ=255)
                                                                       gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAttttt                                                                                                                                        <  4:26749/90‑1 (MQ=255)
                                                                       gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                             >  3:1011304/1‑101 (MQ=255)
                                                                       gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                <  3:911510/146‑1 (MQ=255)
                                                                       gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                >  4:911510/1‑146 (MQ=255)
                                                                       gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaa                                                                             >  1:1251870/1‑149 (MQ=255)
                                                                       gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACGGCATTGATATAAAAAGAAAAAATGATGCCAAAGCGTCATTTTTCCTATCTCAGt                                                                                                                             <  4:1011304/101‑1 (MQ=255)
                                                                        aaaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                          >  2:111424/1‑87 (MQ=38)
                                                                        aaaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  4:517097/149‑1 (MQ=255)
                                                                                tACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTAACAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAAtt                                                                    <  2:1251870/149‑1 (MQ=255)
                                                                                  cAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            >  4:115958/1‑139 (MQ=255)
                                                                                            tGGGCATAATGATGATTGATTTTACTGCATTGATATAAAACGAAAAAATGATGCCAAAGCATCATTTTTc                                                                                                                                       <  3:1195818/70‑1 (MQ=255)
                                                                                            tGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                             <  1:828786/80‑1 (MQ=255)
                                                                                                 aTAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCt                                                                                                                                     >  2:663126/1‑67 (MQ=255)
                                                                                                  tAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  1:887958/1‑52 (MQ=255)
                                                                                                   aaTGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTc                                                                                  <  1:369620/116‑1 (MQ=255)
                                                                                                   aaTGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTc                                                                                  >  2:369620/1‑116 (MQ=255)
                                                                                                    atgatgATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtttt                                                                                                                                         >  2:829311/1‑60 (MQ=255)
                                                                                                    atgatgATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaaa                                                                                           <  3:639485/106‑1 (MQ=255)
                                                                                                     tgatgaTTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcc                                               <  2:1077247/149‑1 (MQ=255)
                                                                                                      gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                           <  3:416182/56‑1 (MQ=25)
                                                                                                      gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                           >  4:416182/1‑56 (MQ=25)
                                                                                                      gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  1:1148314/119‑1 (MQ=255)
                                                                                                       atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  2:965816/118‑1 (MQ=255)
                                                                                                          attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGtt                                                                                                                            >  2:127423/1‑67 (MQ=255)
                                                                                                          attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGtt                                                                                                                            <  1:127423/67‑1 (MQ=255)
                                                                                                          attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGt                                                                                                                         >  2:48574/1‑70 (MQ=255)
                                                                                                          attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGt                                                                                                                         <  1:48574/70‑1 (MQ=255)
                                                                                                          attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATg                                                                                                   >  1:678914/1‑92 (MQ=255)
                                                                                                          attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                          >  3:202915/1‑149 (MQ=255)
                                                                                                           ttgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAAc                                                                                                                 <  4:1202991/77‑1 (MQ=255)
                                                                                                                 ttACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCCGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTa                                   <  4:202915/149‑1 (MQ=255)
                                                                                                                        aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg                                                                                                <  3:1055063/81‑1 (MQ=255)
                                                                                                                        aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                <  2:1200269/97‑1 (MQ=255)
                                                                                                                        aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaaa                                                                       <  1:634487/106‑1 (MQ=255)
                                                                                                                        aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGAAg                                                                                                >  4:1055063/1‑81 (MQ=255)
                                                                                                                          tGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                  <  1:642965/125‑1 (MQ=255)
                                                                                                                            atatAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaaaa                                                                      >  3:602647/1‑103 (MQ=255)
                                                                                                                             tataAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                >  4:320954/1‑92 (MQ=255)
                                                                                                                               taAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttt                                                    <  2:84259/118‑1 (MQ=255)
                                                                                                                                aaaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCa                                                                <  1:1303898/105‑1 (MQ=255)
                                                                                                                                aaaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCa                                                                >  2:1303898/1‑105 (MQ=255)
                                                                                                                                aaaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcc                                               >  3:753028/1‑122 (MQ=255)
                                                                                                                                aaaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcc                                               <  4:753028/122‑1 (MQ=255)
                                                                                                                                 aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTAc                                                                                                            <  1:1237736/60‑1 (MQ=255)
                                                                                                                                   aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCa                                                                >  3:794797/1‑102 (MQ=255)
                                                                                                                                   aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCa                                                                <  4:794797/102‑1 (MQ=255)
                                                                                                                                    aggaaaaaTGATGCTTTGGCATCATTTTTTCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTc                                                                                  <  1:89632/80‑1 (MQ=25)
                                                                                                                                      gaaaaaTGATGCTTTGGCATCATTTTTTCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                  <  2:594692/112‑1 (MQ=38)
                                                                                                                                        aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  4:487532/85‑1 (MQ=255)
                                                                                                                                        aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            >  3:487532/1‑85 (MQ=255)
                                                                                                                                        aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAt                                                      <  3:1083693/107‑1 (MQ=255)
                                                                                                                                        aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTa                                   <  1:833816/126‑1 (MQ=255)
                                                                                                                                         aaaTGATGCTTTGGCATCATTTTTTCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAAtt                                                                    <  4:626293/92‑1 (MQ=38)
                                                                                                                                         aaaTGATGCTTTGGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaaaa                                                                      >  2:457868/1‑90 (MQ=39)
                                                                                                                                         aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                     <  4:373540/123‑1 (MQ=255)
                                                                                                                                         aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCCCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                     >  3:373540/1‑123 (MQ=255)
                                                                                                                                           atgatgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAAt             >  3:1227654/1‑145 (MQ=255)
                                                                                                                                           atgatgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAAt             <  4:1227654/145‑1 (MQ=255)
                                                                                                                                               tgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                  >  4:570703/1‑104 (MQ=255)
                                                                                                                                                gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaa                                                                        <  4:689538/81‑1 (MQ=255)
                                                                                                                                                gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaa                                                                        >  3:689538/1‑81 (MQ=255)
                                                                                                                                                gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGaa         >  1:195330/1‑144 (MQ=255)
                                                                                                                                                  cAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCAtttt  >  2:247620/1‑149 (MQ=255)
                                                                                                                                                  cAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTAGTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAag                          <  1:96807/125‑1 (MQ=255)
                                                                                                                                                    aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATc            >  4:137325/1‑137 (MQ=255)
                                                                                                                                                    aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATc            <  3:137325/137‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTT  >  NC_000913/3792619‑3792913

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