Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 248,769 G→C 8.5% pseudogene (1302/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,7690GC 8.5% 313.6 / 14.2 102pseudogene (1302/1713 nt)lfhApseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (53/40);  new base C (4/5);  total (57/45)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 105 total aligned reads displayed.

ACCACGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATG  >  NC_000913/248620‑248918
                                                                                                                                                     |                                                                                                                                                     
accacGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTg                                                                                                                                                       >  2:1204061/1‑150 (MQ=255)
accacGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTg                                                                                                                                                       >  1:1434178/1‑150 (MQ=255)
   acGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCg                                                                                                                                                    >  2:957336/1‑150 (MQ=255)
   acGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCg                                                                                                                                                    <  1:1412017/150‑1 (MQ=255)
        ggTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGCCCAg                                                                                                                                               >  2:386603/1‑150 (MQ=255)
         gTCAGCAGATTTTCCAGCTCATTTTTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGc                                                                                                                                              >  1:1488689/1‑150 (MQ=255)
          tCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCa                                                                                                                                             >  2:737743/1‑150 (MQ=255)
           cagcagATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCAc                                                                                                                                            >  2:200050/1‑150 (MQ=255)
           cagcagATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCAc                                                                                                                                            >  2:133586/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)
                                           acacCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGaaa                                                                                                            >  1:1463684/1‑150 (MQ=255)
                                             acCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCATAATATGATCTTTCGTCACCGTGCGACTGTCGACCACCACTGTGGCGATGGTGACGATATCGCGCAGGGAAACg                                                                                                          >  2:339286/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                                                                                                   >  1:443271/1‑150 (MQ=255)
                                                    gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCgg                                                                                                   >  2:1260090/1‑150 (MQ=255)
                                                    gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCgg                                                                                                   <  2:1616981/150‑1 (MQ=255)
                                                    gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCgg                                                                                                   >  1:847422/1‑150 (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:957336/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:949308/1‑150 (MQ=255)
                                                                                           gcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAt                                                            <  2:1702496/150‑1 (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)
                                                                                                               gccgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCt                                        >  2:194694/1‑150 (MQ=255)
                                                                                                               gccgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCt                                        >  2:1664583/1‑150 (MQ=255)
                                                                                                                   ccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCg                                    <  2:632001/150‑1 (MQ=255)
                                                                                                                   ccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCg                                    <  2:1495093/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: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)
                                                                                                                                   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)
                                                                                                                                                     |                                                                                                                                                     
ACCACGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATG  >  NC_000913/248620‑248918

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