Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 248,767 C→G 5.1% pseudogene (1304/1713 nt) lfhA ← flagellar system protein, promoterless fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913248,7670CG 5.1% 268.9 / 4.2 78pseudogene (1304/1713 nt)lfhAflagellar system protein, promoterless fragment
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base G (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTT  >  NC_000913/248624‑248915
                                                                                                                                               |                                                                                                                                                    
cGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGa                                                                                                                                                >  1:880995/1‑150 (MQ=255)
  ttGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGAcc                                                                                                                                              <  2:782870/150‑1 (MQ=255)
         gcagATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACgg                                                                                                                                       >  2:156843/1‑150 (MQ=255)
         gcagATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACgg                                                                                                                                       >  2:822999/1‑150 (MQ=255)
           agATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTg                                                                                                                                     <  1:1410547/150‑1 (MQ=255)
             aTTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGAGCTTTCGCCACCGTGCTTGTCGCCAGTAGCACGGTgac                                                                                                                                   <  2:257171/150‑3 (MQ=255)
                  ccAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATgg                                                                                                                              >  1:978667/1‑150 (MQ=255)
                  ccAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGCCACCGTGCTGGTCGCCAGTAGCACGGTGACGaaag                                                                                                                              <  2:132354/150‑4 (MQ=38)
                   cAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGt                                                                                                                             >  1:770304/1‑150 (MQ=255)
                     gCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGa                                                                                                                           >  1:85295/1‑150 (MQ=255)
                       tCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACg                                                                                                                         <  2:1324374/150‑1 (MQ=255)
                        cATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGAc                                                                                                                                               <  2:119643/127‑1 (MQ=255)
                        cATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGAc                                                                                                                                               >  1:119643/1‑127 (MQ=255)
                        cATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGa                                                                                                                        >  2:571994/1‑150 (MQ=255)
                        cATTATTCAGCGTATACACCATCAGCTCCTGCTTGAGAACGAACGGATCGGTAATGCTGCGCCGCAATGCCAGGCGCAAATCGGACGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGa                                                                                                                        >  2:484851/1‑150 (MQ=255)
                             ttCAGCGTATACACCGTAAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATcg                                                                                                                   >  1:450353/1‑150 (MQ=255)
                               cAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATcgcg                                                                                                                 <  1:620531/150‑1 (MQ=255)
                                 gCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCa                                                                                                               <  1:450355/150‑1 (MQ=255)
                                     ataCACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGa                                                                                                           >  2:540094/1‑150 (MQ=255)
                                       acacCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGaaa                                                                                                         >  2:924833/1‑150 (MQ=255)
                                         acCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACg                                                                                                       <  2:888347/150‑1 (MQ=255)
                                                gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCgg                                                                                                <  1:1032232/150‑1 (MQ=255)
                                                gCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATACCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCgg                                                                                                >  1:627972/1‑150 (MQ=255)
                                                   ccTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTca                                                                                             <  2:746083/150‑1 (MQ=255)
                                                     tGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTcagc                                                                                           <  2:1176214/150‑1 (MQ=255)
                                                         tGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAgcg                                                                                       <  1:156843/150‑1 (MQ=255)
                                                           cgaacgaaCGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGAc                                                                                                                          >  1:1204772/1‑113 (MQ=255)
                                                           cgaacgaaCGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGAc                                                                                                                          <  2:1204772/113‑1 (MQ=255)
                                                           cgaacgaaCGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGcc                                                                                     >  1:1386249/1‑150 (MQ=255)
                                                               cgaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGt                                                                                 <  2:990371/150‑1 (MQ=255)
                                                               cgaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGt                                                                                 >  2:1354434/1‑150 (MQ=255)
                                                               cgaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGt                                                                                 <  2:1408172/150‑1 (MQ=255)
                                                                gaacgGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTa                                                                                <  1:430095/150‑1 (MQ=255)
                                                                        gggTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCa                                                                        <  2:56418/150‑1 (MQ=255)
                                                                        gggTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCa                                                                        <  1:924833/150‑1 (MQ=255)
                                                                        gggTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCa                                                                        <  2:79989/150‑1 (MQ=255)
                                                                        gggTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCa                                                                        <  2:978667/150‑1 (MQ=255)
                                                                              tgctgcGCCGCAATGCCAGGCGCACATCGGTCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACt                                                                  >  1:484871/1‑150 (MQ=255)
                                                                               gctgcGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTg                                                                 >  2:650326/1‑150 (MQ=255)
                                                                                 tgcGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGc                                                               <  2:768311/150‑1 (MQ=255)
                                                                                   cgccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGTCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTg                                                             <  2:509427/150‑1 (MQ=255)
                                                                                   cgccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTg                                                             >  2:1030069/1‑150 (MQ=255)
                                                                                    gccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGt                                                            <  1:540094/150‑1 (MQ=255)
                                                                                    gccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGt                                                            >  1:855683/1‑150 (MQ=255)
                                                                                     ccgcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTa                                                           >  1:1238040/1‑150 (MQ=255)
                                                                                       gcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAt                                                         >  1:234767/1‑150 (MQ=255)
                                                                                       gcAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAt                                                         >  1:1051295/1‑150 (MQ=255)
                                                                                         aaTGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTg                                                       <  2:1359094/150‑1 (MQ=255)
                                                                                            gCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAg                                                                       >  1:432071/1‑131 (MQ=255)
                                                                                            gCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAg                                                                       <  2:432071/131‑1 (MQ=255)
                                                                                            gCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAgc                                                    <  1:91640/150‑1 (MQ=255)
                                                                                                   gcACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAt                                                         <  1:36818/138‑1 (MQ=255)
                                                                                                   gcACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAAt                                                         >  2:36818/1‑138 (MQ=255)
                                                                                                     acaTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTc                                           >  1:1412894/1‑150 (MQ=255)
                                                                                                        tCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaaa                                        >  1:290092/1‑150 (MQ=255)
                                                                                                        tCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaaa                                        <  1:660859/150‑1 (MQ=255)
                                                                                                        tCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaaa                                        >  2:635342/1‑150 (MQ=255)
                                                                                                          ggCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCGGCGCTCaatta                                      >  2:712077/1‑147 (MQ=255)
                                                                                                             cgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCGGCGCTCaattacag                                   <  1:712077/150‑7 (MQ=255)
                                                                                                              gccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTTTAATTGAGCGCGGCGCTCaattacagc                                  <  2:450353/150‑8 (MQ=255)
                                                                                                               ccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCg                                 <  2:1238040/150‑1 (MQ=255)
                                                                                                               ccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCg                                 <  2:122645/150‑1 (MQ=255)
                                                                                                                 agcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGcc                               <  2:1051295/150‑1 (MQ=255)
                                                                                                                  gcagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTAGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCa                              >  2:207386/1‑150 (MQ=255)
                                                                                                                        tatGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgcg                        >  2:1370380/1‑150 (MQ=255)
                                                                                                                         atGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTTGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCgg                       >  2:551809/1‑150 (MQ=255)
                                                                                                                         atGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCgg                       <  1:1370380/150‑1 (MQ=255)
                                                                                                                           gATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTg                     >  1:567979/1‑150 (MQ=255)
                                                                                                                           gATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTg                     >  1:647226/1‑150 (MQ=255)
                                                                                                                               tttCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGt                 >  1:1380021/1‑150 (MQ=255)
                                                                                                                                   gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac             <  1:484851/150‑1 (MQ=255)
                                                                                                                                   gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac             <  2:484871/150‑1 (MQ=255)
                                                                                                                                    ccACCGTGCTGGTCGCCAGTAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCCCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgacg            >  2:1009613/2‑150 (MQ=255)
                                                                                                                                    ccACCGTGCTGGTCGCCAGTAGCACGGTGACGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgacg            >  1:1074693/2‑150 (MQ=255)
                                                                                                                                     cACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTTGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgacga           <  1:551809/150‑1 (MQ=255)
                                                                                                                                     cACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgacga           <  1:571994/150‑1 (MQ=255)
                                                                                                                                            cTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACgg    >  1:1315858/1‑150 (MQ=255)
                                                                                                                                              aCTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTTTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGtt  <  2:838291/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CGTTGGTCAGCAGATTTTCCAGCTCATTATTCAGCGTATACACCGTCAGCTCCTGCTTGCGAACGAACGGATGGGTAATGCTGCGCCGCAATGCCAGGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTT  >  NC_000913/248624‑248915

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