Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,792,766 A→T 8.2% intergenic (‑215/+60) 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,7660AT 8.2% 601.5 / 32.4 171intergenic (‑215/+60)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 = 6.10e-01

Only 100 of 184 total aligned reads displayed.

GCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTT  >  NC_000913/3792619‑3792913
                                                                                                                                                   |                                                                                                                                                   
          ttCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                          <  1:828530/149‑1 (MQ=255)
          ttCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACAACTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                          <  1:1455494/149‑1 (MQ=255)
               cAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCt                                                                                                                                     >  4:837571/1‑149 (MQ=255)
                    ttCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTAtctc                                                                                                                                >  3:757743/1‑149 (MQ=255)
                        ggATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtttt                                                                                                                                         >  4:690072/1‑136 (MQ=255)
                         gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  3:1070572/1‑125 (MQ=255)
                         gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  3:1437046/1‑125 (MQ=255)
                         gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  4:295250/125‑1 (MQ=255)
                         gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTg                                                                                                                           >  1:1025933/1‑149 (MQ=255)
                               taaCGTGTAACTCGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtttt                                                                                                                                         <  2:1330391/129‑1 (MQ=255)
                               taaCGTGTAACTCGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtttt                                                                                                                                         >  1:1330391/1‑129 (MQ=255)
                                    tgtAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  4:714861/1‑114 (MQ=255)
                                    tgtAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  3:714861/114‑1 (MQ=255)
                                                  tctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatcat                                                                                                                                            >  2:308424/1‑107 (MQ=255)
                                                  tctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATgg                                                                                                  >  1:1077247/1‑149 (MQ=255)
                                                         acCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcat                                                                                                                                               <  1:767206/97‑1 (MQ=255)
                                                         acCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcat                                                                                                                                               >  2:767206/1‑97 (MQ=255)
                                                          ccTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAAc                                                                                          <  3:318209/149‑1 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatca                                                                                                                                             >  3:271174/1‑97 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAttt                                                                                                                                          >  2:1111867/1‑100 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTg                                                                                                                           <  4:1403778/115‑1 (MQ=255)
                                                           cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAAc                                                                                          >  4:396438/1‑148 (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)
                                                                     tagAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGt                                                                                                                         <  2:567825/107‑1 (MQ=255)
                                                                       gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAttttt                                                                                                                                        <  4:26749/90‑1 (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)
                                                                        aaaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  4:517097/149‑1 (MQ=255)
                                                                        aaaaaCCGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTc                                                                                  <  3:809095/143‑1 (MQ=255)
                                                                                tACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAgg                                                                                               <  3:1025968/122‑1 (MQ=255)
                                                                                tACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTAACAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAAtt                                                                    <  2:1251870/149‑1 (MQ=255)
                                                                                  cAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  3:115958/139‑1 (MQ=255)
                                                                                            tGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAACTGATGCCAAAGCATCATTTTTc                                                                                                                                       >  4:1195818/1‑70 (MQ=255)
                                                                                            tGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                             >  2:828786/1‑80 (MQ=255)
                                                                                                  tAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  3:1344225/1‑52 (MQ=255)
                                                                                                  tAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  2:887958/52‑1 (MQ=255)
                                                                                                  tAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  1:887958/1‑52 (MQ=255)
                                                                                                   aaTGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTc                                                                                  <  1:369620/116‑1 (MQ=255)
                                                                                                    atgatgATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                          >  1:560579/1‑59 (MQ=25)
                                                                                                    atgatgATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtttt                                                                                                                                         <  1:829311/60‑1 (MQ=255)
                                                                                                    atgatgATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaaa                                                                                           <  3:639485/106‑1 (MQ=255)
                                                                                                     tgatgaTTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcc                                               <  2:1077247/149‑1 (MQ=255)
                                                                                                      gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                           <  3:416182/56‑1 (MQ=25)
                                                                                                      gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                <  4:1064247/115‑1 (MQ=255)
                                                                                                      gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                >  3:1064247/1‑115 (MQ=255)
                                                                                                      gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  1:1148314/119‑1 (MQ=255)
                                                                                                      gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            >  2:1148314/1‑119 (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)
                                                                                                          attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATg                                                                                                   <  2:678914/92‑1 (MQ=255)
                                                                                                          attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                          >  3:202915/1‑149 (MQ=255)
                                                                                                           ttgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAAc                                                                                                                 >  3:1202991/1‑77 (MQ=255)
                                                                                                           ttgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTAtt                                                                                     <  2:608423/105‑1 (MQ=255)
                                                                                                           ttgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTAtt                                                                                     >  1:608423/1‑105 (MQ=255)
                                                                                                           ttgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGc                                                        <  4:177147/134‑1 (MQ=255)
                                                                                                                 ttACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCCGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTa                                   <  4:202915/149‑1 (MQ=255)
                                                                                                                        aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaaa                                                                       <  1:634487/106‑1 (MQ=255)
                                                                                                                        aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGAAg                                                                                                >  4:1055063/1‑81 (MQ=255)
                                                                                                                        aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGGTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaaa                                                                       >  2:634487/1‑106 (MQ=255)
                                                                                                                          tGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                  >  2:642965/1‑125 (MQ=255)
                                                                                                                          tGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                  <  1:642965/125‑1 (MQ=255)
                                                                                                                             tataAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                >  4:320954/1‑92 (MQ=255)
                                                                                                                              ataAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgt                                       >  2:392519/1‑132 (MQ=255)
                                                                                                                               taAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttt                                                    >  1:84259/1‑118 (MQ=255)
                                                                                                                                aaaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCa                                                                <  1:1303898/105‑1 (MQ=255)
                                                                                                                                aaaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCa                                                                >  2:1303898/1‑105 (MQ=255)
                                                                                                                                aaaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcc                                               <  4:753028/122‑1 (MQ=255)
                                                                                                                                aaaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgt                                       <  3:250343/130‑1 (MQ=255)
                                                                                                                                 aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTAc                                                                                                            >  2:1237736/1‑60 (MQ=255)
                                                                                                                                 aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTAc                                                                                                            <  1:1237736/60‑1 (MQ=255)
                                                                                                                                 aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTc                                                                                  >  3:751976/1‑86 (MQ=255)
                                                                                                                                   aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                                 <  4:485501/85‑1 (MQ=255)
                                                                                                                                   aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCa                                                                <  4:794797/102‑1 (MQ=255)
                                                                                                                                    aggaaaaaTGATGCTTTGGCATCATTTTTTCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTc                                                                                  <  1:89632/80‑1 (MQ=25)
                                                                                                                                    aggaaaaaTGATGCTTTGGCATCATTTTTTCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTc                                                                                  >  2:89632/4‑83 (MQ=25)
                                                                                                                                        aaaaTGATGCCAAAGTATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAt                                                      >  4:1083693/1‑107 (MQ=255)
                                                                                                                                        aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            >  3:487532/1‑85 (MQ=255)
                                                                                                                                        aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  4:487532/85‑1 (MQ=255)
                                                                                                                                        aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTa                                   >  2:833816/1‑126 (MQ=255)
                                                                                                                                         aaaTGATGCTTTGGCATCATTTTTTCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAAtt                                                                    <  4:626293/92‑1 (MQ=38)
                                                                                                                                         aaaTGATGCTTTGGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaaaa                                                                      >  2:457868/1‑90 (MQ=39)
                                                                                                                                         aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCCCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                     >  3:373540/1‑123 (MQ=255)
                                                                                                                                           atgatgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAAt             <  4:1227654/145‑1 (MQ=255)
                                                                                                                                               tgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                  >  4:570703/1‑104 (MQ=255)
                                                                                                                                               tgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                  <  3:570703/104‑1 (MQ=255)
                                                                                                                                                gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaa                                                                        >  3:689538/1‑81 (MQ=255)
                                                                                                                                                gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGaa         >  1:195330/1‑144 (MQ=255)
                                                                                                                                                gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGaa         <  2:195330/144‑1 (MQ=255)
                                                                                                                                                 ccAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTAc                                  >  3:495629/1‑118 (MQ=255)
                                                                                                                                                  cAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAag                          >  2:96807/1‑125 (MQ=255)
                                                                                                                                                  cAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTAGTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAag                          <  1:96807/125‑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 < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: