Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 248,772 G→C 6.6% pseudogene (1299/1713 nt) lfhA ← pseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913248,7720GC 6.6% 341.4 / 10.1 103pseudogene (1299/1713 nt)lfhApseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (55/41);  new base C (3/4);  total (58/45)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

ACGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAA  >  NC_000913/248623‑248921
                                                                                                                                                     |                                                                                                                                                     
acGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCg                                                                                                                                                       >  2:957336/1‑150 (MQ=255)
acGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCg                                                                                                                                                       <  1:1412017/150‑1 (MQ=255)
     ggTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGCCCAg                                                                                                                                                  >  2:386603/1‑150 (MQ=255)
       tCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCa                                                                                                                                                >  2:737743/1‑150 (MQ=255)
        cagcagATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCAc                                                                                                                                               >  2:133586/1‑150 (MQ=255)
        cagcagATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCAc                                                                                                                                               >  2:200050/1‑150 (MQ=255)
        cagcagATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGAGCTTTCGTCACCGTGCTACTGGCGACCAGCAc                                                                                                                                               >  2:1609881/1‑150 (MQ=255)
          gcagATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACgg                                                                                                                                             >  2:780510/1‑150 (MQ=255)
              aTTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGCGACCGTGCTGGTCGCCAGTAGCACGGTgac                                                                                                                                         <  1:1501181/150‑3 (MQ=255)
              aTTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGCCACCGTGCTGGTCGCCAGTAGCACGGTgac                                                                                                                                         <  1:1229812/150‑3 (MQ=255)
               ttttCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCg                                                                                                                                        <  2:12506/150‑1 (MQ=255)
                  tCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCTGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATg                                                                                                                                     <  2:559550/150‑1 (MQ=255)
                   ccAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATgg                                                                                                                                    >  1:1182655/1‑150 (MQ=255)
                   ccAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGCCACCGTGCTGGTCGCCAGTAGCACGGTGACGaaag                                                                                                                                    <  1:524783/150‑4 (MQ=38)
                   ccAGCTCATATTTCACTGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGATTATGATCTTTCGCCACCGGGCTGGTCGCCAGTAGCACGGTGACGaaag                                                                                                                                    <  2:186404/150‑4 (MQ=37)
                      gCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGAGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGa                                                                                                                                 >  2:514543/1‑150 (MQ=255)
                       cTCATTATTCATCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGCATGGGTAATGCTGCGCCGCAAGGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGGCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGAc                                                                                                                                <  2:1653798/150‑1 (MQ=255)
                       cTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAACGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGAc                                                                                                                                >  1:635703/1‑150 (MQ=255)
                         cATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGa                                                                                                                              <  2:640745/150‑1 (MQ=255)
                            tattCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGCTCTTTCGTCACAGTGCTACTCGCGACCAGCACGGTCGCGATGGTGACGatat                                                                                                                           >  1:68534/1‑150 (MQ=255)
                             attCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATc                                                                                                                          >  2:690697/1‑150 (MQ=255)
                                     tataCACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATCTCGCGCAggg                                                                                                                  >  1:1525562/1‑150 (MQ=255)
                                        acacCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGaaa                                                                                                               >  1:1463684/1‑150 (MQ=255)
                                          acCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGCTGGTGACGATATCGCGCAGGGAAACg                                                                                                             >  1:1011656/1‑150 (MQ=255)
                                          acCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACg                                                                                                             <  2:561195/150‑1 (MQ=255)
                                          acCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAAAGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTTGGGATGGTGACGATATCGCGCAGGGAAACg                                                                                                             >  2:334205/1‑150 (MQ=255)
                                             gTCAGCTCCTGCGTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCt                                                                                                          <  1:1476373/150‑1 (MQ=255)
                                               cAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTc                                                                                                        <  1:1592216/150‑1 (MQ=255)
                                                aGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCg                                                                                                       >  1:1146087/1‑150 (MQ=255)
                                                 gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCgg                                                                                                      >  2:1260090/1‑150 (MQ=255)
                                                 gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCgg                                                                                                      >  1:847422/1‑150 (MQ=255)
                                                 gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCgg                                                                                                      >  1:443271/1‑150 (MQ=255)
                                                 gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCgg                                                                                                      <  2:1616981/150‑1 (MQ=255)
                                                       gCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGAGATCGCGCAGGGAACCGCCTTCGGTcagcc                                                                                                >  2:285361/1‑149 (MQ=255)
                                                        cTTGCGAACGAACGGAGGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTATCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTcagcag                                                                                               <  2:1724718/150‑1 (MQ=255)
                                                          tGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAgcg                                                                                             >  2:6852/1‑150 (MQ=255)
                                                               acgaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGGGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCgg                                                                                        <  1:1204061/150‑1 (MQ=255)
                                                                    cgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTacac                                                                                   >  1:301813/1‑150 (MQ=255)
                                                                     ggATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACt                                                                                  <  1:554839/150‑1 (MQ=255)
                                                                      gATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACtt                                                                                 >  2:333589/1‑150 (MQ=255)
                                                                         gggTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCa                                                                              <  1:924478/150‑1 (MQ=255)
                                                                           gTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGc                                                                            <  1:188197/150‑1 (MQ=255)
                                                                              aTGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAAc                                                                         >  1:1528585/1‑150 (MQ=255)
                                                                                   gcgcCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCt                                                                    >  1:616337/1‑150 (MQ=255)
                                                                                    cgccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGGCACCGTGCTACTGGGGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTg                                                                   <  1:243577/150‑1 (MQ=255)
                                                                                    cgccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGAGCTTTCGGCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTg                                                                   <  1:350659/150‑1 (MQ=255)
                                                                                     gccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGt                                                                  >  1:1412809/1‑150 (MQ=255)
                                                                                      ccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTa                                                                 >  2:1422675/1‑150 (MQ=255)
                                                                                        gcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAt                                                               <  2:1702496/150‑1 (MQ=255)
                                                                                        gcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAt                                                               >  2:949308/1‑150 (MQ=255)
                                                                                        gcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAt                                                               >  1:23759/1‑150 (MQ=255)
                                                                                           cattgcggCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGa                                                            >  1:587227/7‑150 (MQ=255)
                                                                                            tGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAg                                                           <  2:1146087/150‑1 (MQ=255)
                                                                                              ccAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCCGCCCGGTGGCGATGGTGACGATATCGCGCAGGGACACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAgcg                                                         >  1:883593/1‑150 (MQ=255)
                                                                                              ccAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAgcg                                                         <  1:133586/150‑1 (MQ=255)
                                                                                               cAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAgcgc                                                        >  1:526064/1‑150 (MQ=255)
                                                                                                 ggCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGcgccg                                                      >  1:59774/1‑150 (MQ=255)
                                                                                                   cgcACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCcgcg                                                    <  1:1285532/150‑1 (MQ=255)
                                                                                                    gcACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCgcgc                                                   >  2:1165384/1‑150 (MQ=255)
                                                                                                        aTCGGCCTCCAGCAGAATATGATCTTTCGGCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGAGATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCGGTAATTGAGCGCCGCGCTCaa                                               <  2:1094167/150‑1 (MQ=255)
                                                                                                        aTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaa                                               >  1:1333690/1‑150 (MQ=255)
                                                                                                          cGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAAt                                             >  2:325712/1‑150 (MQ=255)
                                                                                                            gccgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCGGCGCTCaattac                                           <  2:301813/150‑5 (MQ=255)
                                                                                                            gccgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCt                                           >  2:194694/1‑150 (MQ=255)
                                                                                                            gccgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCt                                           >  2:1664583/1‑150 (MQ=255)
                                                                                                                ccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCg                                       <  2:632001/150‑1 (MQ=255)
                                                                                                                 cagcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGc                                      <  2:1409672/150‑1 (MQ=255)
                                                                                                                  agcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGcc                                     <  2:1027216/150‑1 (MQ=255)
                                                                                                                   gcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCa                                    <  2:1359149/150‑1 (MQ=255)
                                                                                                                     agAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGa                                  <  2:1300846/150‑1 (MQ=255)
                                                                                                                     agAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGa                                  <  2:1397995/150‑1 (MQ=255)
                                                                                                                     agAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGa                                  >  1:1406513/1‑150 (MQ=255)
                                                                                                                      gAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCTGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAc                                 <  1:962551/150‑1 (MQ=255)
                                                                                                                        atatGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgc                               >  2:657868/1‑150 (MQ=255)
                                                                                                                        atatGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgc                               >  2:28239/1‑150 (MQ=255)
                                                                                                                           tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGCCGCGGt                            >  2:430721/1‑150 (MQ=255)
                                                                                                                           tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGt                            <  1:1664262/150‑1 (MQ=255)
                                                                                                                            gATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTg                           <  1:200050/150‑1 (MQ=255)
                                                                                                                            gATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTg                           <  1:1609881/150‑1 (MQ=255)
                                                                                                                            gATCTTGCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTGCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTg                           <  1:418231/150‑1 (MQ=255)
                                                                                                                             ctCTTTCGTCTCCGTCCTAGTCGCGACCAGCACGTTGTCGATGGTGACGATATCGCGCTGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCTCGCTCAAATCTTCCGCCAGACGCGTTGc                          <  2:277275/149‑1 (MQ=255)
                                                                                                                             aTCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGc                          <  1:386603/150‑1 (MQ=255)
                                                                                                                               ccatCGCCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCCGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCg                        >  1:671903/4‑150 (MQ=255)
                                                                                                                               cTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCg                        >  1:606860/1‑150 (MQ=255)
                                                                                                                                tttCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGt                       >  1:1002231/1‑150 (MQ=255)
                                                                                                                                 ttCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTc                      <  1:325712/150‑1 (MQ=255)
                                                                                                                                 atCGCCACCGTGCTGGTCGCCAGTAGCACGGTGACGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGCGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTc                      >  1:938139/2‑150 (MQ=38)
                                                                                                                                  tCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcg                     <  2:616337/150‑1 (MQ=255)
                                                                                                                                   cgccACCGTGCTGGTCGCCAGTAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcga                    >  2:74654/4‑150 (MQ=255)
                                                                                                                                   cgccACCGTGCTGGTCGCCAGTAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTACGCCAGACGCGGTGCCGTcga                    >  2:1266098/4‑150 (MQ=255)
                                                                                                                                   cGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgc                    >  2:1268916/1‑149 (MQ=255)
                                                                                                                                       aCCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGAt                >  1:1184014/1‑150 (MQ=255)
                                                                                                                                         cGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATaa              <  2:1434178/150‑1 (MQ=255)
                                                                                                                                           tGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAAc            >  2:1508476/1‑150 (MQ=255)
                                                                                                                                            gCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCCGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACg           >  1:20440/1‑150 (MQ=255)
                                                                                                                                            gCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACg           <  1:690697/150‑1 (MQ=255)
                                                                                                                                               aCTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCCGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGtt        >  1:1135754/1‑150 (MQ=255)
                                                                                                                                                  ggCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATg     >  1:1425335/1‑150 (MQ=255)
                                                                                                                                                   gCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGc    <  2:1221097/150‑1 (MQ=255)
                                                                                                                                                    cGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCa   >  1:936781/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAA  >  NC_000913/248623‑248921

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