Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 3,230,291 (T)6→7 intergenic (‑112/+961) rfaY ← / ← rfaI lipopolysaccharide core biosynthesis protein/UDP‑D‑galactose:(glucosyl)lipopolysaccharide‑ alpha‑1,3‑D‑galactosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123063,230,2851.T96.9% 1882.1 / 40.6 481intergenic (‑106/+967)rfaY/rfaIlipopolysaccharide core biosynthesis protein/UDP‑D‑galactose:(glucosyl)lipopolysaccharide‑ alpha‑1,3‑D‑galactosyltransferase
Reads supporting (aligned to +/- strand):  ref base . (6/9);  new base T (246/220);  total (252/229)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 528 total aligned reads displayed.

TCGGTAAAAACAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAATGAATAAG  >  AP012306/3230136‑3230429
                                                                                                                                                      |                                                                                                                                                
tCGGTAAAAACAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATc                                                                                                                                                   <  1:2492756/150‑1 (MQ=255)
 cGGTAAAAACAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                                                                                                                                 <  1:3556121/150‑1 (MQ=255)
 cGGTAAAAACAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                                                                                                                                 <  1:4047255/150‑1 (MQ=255)
    tAAAAACAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                                                                                                                              <  2:3371153/150‑1 (MQ=255)
    tAAAAACAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                                                                                                                              <  1:4963713/150‑1 (MQ=255)
     aaaaaCAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                                                                                                                             >  2:2048766/1‑150 (MQ=255)
     aaaaaCAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                                                                                                                             >  1:2036902/1‑150 (MQ=255)
       aaaCAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttt                                                                                                                                           <  1:2908345/150‑2 (MQ=255)
       aaaCAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttt                                                                                                                                           <  1:4461965/150‑2 (MQ=255)
        aacaacCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttta                                                                                                                                          >  1:1798534/1‑148 (MQ=255)
            acCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAAtt                                                                                                                                       <  1:706111/150‑1 (MQ=255)
                aGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCaat                                                                                                                                   <  1:544590/150‑1 (MQ=255)
                aGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCaat                                                                                                                                   >  2:459855/1‑150 (MQ=255)
                   cTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCaataat                                                                                                                                >  1:273562/1‑150 (MQ=255)
                    tttGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAAtt                                                                                                                               >  2:1445101/1‑150 (MQ=255)
                      tGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGa                                                                                                                             >  1:3844807/1‑150 (MQ=255)
                       gATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGAt                                                                                                                            <  2:2669524/150‑1 (MQ=255)
                        aTCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGAtt                                                                                                                           <  2:1205119/150‑1 (MQ=255)
                        aTCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGAtt                                                                                                                           <  1:1537351/150‑1 (MQ=255)
                        aTCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGAtt                                                                                                                           <  2:319097/150‑1 (MQ=255)
                         tCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGAttt                                                                                                                          >  2:3838091/1‑150 (MQ=255)
                           ttGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAGTAATTGATTTCg                                                                                                                        <  1:1837600/150‑1 (MQ=255)
                            tGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGc                                                                                                                       <  2:880600/150‑1 (MQ=255)
                            tGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGc                                                                                                                       <  1:3397122/150‑1 (MQ=255)
                             gCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCa                                                                                                                      <  2:248/150‑1 (MQ=255)
                             gCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCa                                                                                                                      <  2:4038087/150‑1 (MQ=255)
                               tcttctGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCAtc                                                                                                                    >  2:1787946/1‑150 (MQ=255)
                               tcttctGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCAtc                                                                                                                    <  2:2147213/150‑1 (MQ=255)
                                 ttctGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCAtctc                                                                                                                  >  1:872561/1‑150 (MQ=255)
                                   ctGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGt                                                                                                                >  1:2242127/1‑150 (MQ=255)
                                      aaTCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga                                                                                                             >  1:3770338/1‑150 (MQ=255)
                                       aTCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGgag                                                                                                            >  1:1015892/1‑150 (MQ=255)
                                        tCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGgaga                                                                                                           >  1:3968681/1‑150 (MQ=255)
                                        tCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGgaga                                                                                                           >  2:286969/1‑150 (MQ=255)
                                            aataaACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa                                                                                                       >  1:2535547/1‑150 (MQ=255)
                                            aataaACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa                                                                                                       >  1:3183045/1‑150 (MQ=255)
                                             ataaACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt                                                                                                      >  2:360361/1‑150 (MQ=255)
                                              taaACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc                                                                                                     >  1:257863/1‑150 (MQ=255)
                                               aaaCCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCt                                                                                                    >  1:1715725/1‑150 (MQ=255)
                                                  ccAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTc                                                                                                 >  2:561553/1‑150 (MQ=255)
                                                    aGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCa                                                                                               >  1:2390406/1‑150 (MQ=255)
                                                    aGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCa                                                                                               <  2:1096391/150‑1 (MQ=255)
                                                      ttAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAgg                                                                                             >  1:1607631/1‑150 (MQ=255)
                                                      ttAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAgg                                                                                             >  1:4794346/1‑150 (MQ=255)
                                                        aaaTGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAgggg                                                                                           >  2:1163877/1‑150 (MQ=255)
                                                                  tACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTa                                                                                 >  1:1423003/1‑150 (MQ=255)
                                                                     ggTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAg                                                                              <  2:1387548/150‑1 (MQ=255)
                                                                        aaTATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTa                                                                           <  1:4126817/150‑1 (MQ=255)
                                                                          tatTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTAtt                                                                         >  1:4894010/1‑150 (MQ=255)
                                                                          tatTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTAtt                                                                         >  2:1722317/1‑150 (MQ=255)
                                                                           atTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTAttt                                                                        >  2:3266090/1‑150 (MQ=255)
                                                                            ttttCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTAtttt                                                                       <  1:422888/150‑1 (MQ=255)
                                                                             tttCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTta                                                                      <  2:2211212/150‑1 (MQ=255)
                                                                             tttCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTta                                                                      >  2:4131054/1‑150 (MQ=255)
                                                                             tttCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTta                                                                      >  2:529106/1‑150 (MQ=255)
                                                                                  gCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTa                                                                 <  2:4869930/150‑1 (MQ=255)
                                                                                  gCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTa                                                                 <  2:3444618/150‑1 (MQ=255)
                                                                                  gCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTa                                                                 <  2:2733586/150‑1 (MQ=255)
                                                                                   cAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTAt                                                                >  2:312791/1‑150 (MQ=255)
                                                                                   cAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTAt                                                                >  1:305276/1‑150 (MQ=255)
                                                                                    aaaGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTAtt                                                               >  1:4541775/1‑150 (MQ=255)
                                                                                     aaGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTAttt                                                              <  1:4195067/150‑1 (MQ=255)
                                                                                         tAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTCGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTAcc                                                          <  1:3373628/150‑1 (MQ=255)
                                                                                          aaCAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCg                                                         >  2:71314/1‑150 (MQ=255)
                                                                                           aCAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGa                                                        >  2:1579546/1‑150 (MQ=255)
                                                                                             aGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATg                                                      >  2:1497326/1‑150 (MQ=255)
                                                                                             aGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATg                                                      >  2:4259245/1‑150 (MQ=255)
                                                                                              gACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATgg                                                     >  1:3218090/1‑150 (MQ=255)
                                                                                              gACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATgg                                                     >  2:2353481/1‑150 (MQ=255)
                                                                                               aCTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGa                                                    >  2:4738399/1‑150 (MQ=255)
                                                                                                 tCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATa                                                  <  1:2173883/150‑1 (MQ=255)
                                                                                                   cTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATaaa                                                >  2:3693917/1‑150 (MQ=255)
                                                                                                   cTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATaaa                                                <  1:3300692/150‑1 (MQ=255)
                                                                                                     gCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAAtt                                              <  2:3532601/150‑1 (MQ=255)
                                                                                                             tCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCAttt                                      >  1:1307540/1‑150 (MQ=255)
                                                                                                              ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTtg                                     <  1:4708481/150‑1 (MQ=255)
                                                                                                               cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTtgt                                    >  1:1685681/1‑150 (MQ=255)
                                                                                                                 gagaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGc                                  <  2:860626/150‑1 (MQ=255)
                                                                                                                     tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACg                              >  1:4593828/1‑150 (MQ=255)
                                                                                                                      ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACgg                             <  1:1931587/150‑1 (MQ=255)
                                                                                                                       taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGt                            >  1:4307764/1‑150 (MQ=255)
                                                                                                                       taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGt                            <  2:3954395/150‑1 (MQ=255)
                                                                                                                          tgatgTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTa                         >  1:2996388/1‑150 (MQ=255)
                                                                                                                                   aCTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCt                <  1:986438/150‑1 (MQ=255)
                                                                                                                                    cTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTg               >  2:478995/1‑150 (MQ=255)
                                                                                                                                     tAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGt              >  1:761625/1‑150 (MQ=255)
                                                                                                                                        aaGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGtata           >  2:2834755/1‑150 (MQ=255)
                                                                                                                                         aGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATaa          >  2:776881/1‑150 (MQ=255)
                                                                                                                                          gATGTTTATATCTTTTTTTAAATTTAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAAt         <  2:3047047/150‑1 (MQ=255)
                                                                                                                                          gATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAAt         >  1:3904489/1‑150 (MQ=255)
                                                                                                                                          gATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAAt         <  1:1100070/150‑1 (MQ=255)
                                                                                                                                          gATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTATAGTTGCACCTGTATAAt         >  2:4386883/1‑150 (MQ=255)
                                                                                                                                            tGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAATGa       >  2:4675520/1‑150 (MQ=255)
                                                                                                                                            tGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAATGa       <  1:36538/150‑1 (MQ=255)
                                                                                                                                             gTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAATGaa      <  2:2556400/150‑1 (MQ=255)
                                                                                                                                              tttATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAATGAAt     >  2:3650593/1‑150 (MQ=255)
                                                                                                                                              tttATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAATGAAt     >  2:2997522/1‑150 (MQ=255)
                                                                                                                                              tttATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAATGAAt     >  2:1924755/1‑150 (MQ=255)
                                                                                                                                                 atatCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAATGAATaag  >  2:262321/1‑150 (MQ=255)
                                                                                                                                                 atatCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAATGAATaag  >  1:2084181/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
TCGGTAAAAACAACCAAGTCTTTGATCTTGCTCTTCTGAATCATAATAAACCAGTTAAATGTTATTTACGGTAATATTTTCGGCAAAGATAACAGACTCCTGCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTCCAGGGGGAATTTTCTAAAGCTATTTTATAGTATTTTACCGATGGATAAATTGCCCATTTGTGCCACGGTTTAGTTGCACCTGTATAATGAATAAG  >  AP012306/3230136‑3230429

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: