Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,141,148 C→T 8.7% intergenic (+162/+34) flgL → / ← rne flagellar hook‑filament junction protein 2/ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,141,1480CT 8.7% 271.2 / 10.8 93intergenic (+162/+34)flgL/rneflagellar hook‑filament junction protein 2/ribonuclease E
Reads supporting (aligned to +/- strand):  ref base C (41/43);  new base T (4/4);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CATATCATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGA  >  NC_000913/1141000‑1141296
                                                                                                                                                    |                                                                                                                                                    
cATATCATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCa                                                                                                                                                     >  1:348388/1‑150 (MQ=255)
cATATCATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCa                                                                                                                                                     >  2:1038938/1‑150 (MQ=255)
     cATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTAc                                                                                                                                                >  2:615672/1‑150 (MQ=255)
     cATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTAc                                                                                                                                                >  2:522013/1‑150 (MQ=255)
         aaaCTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAgg                                                                                                                                            >  2:873889/1‑150 (MQ=255)
           aCTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGc                                                                                                                                          <  1:792827/150‑1 (MQ=255)
            cTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCt                                                                                                                                         <  2:1126990/150‑1 (MQ=255)
             tGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGTAACTGCCAGGGCtt                                                                                                                                        <  1:683641/150‑1 (MQ=255)
              gggTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTg                                                                                                                                       <  1:511029/150‑1 (MQ=255)
               ggTATGTTTTGTCTACCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGa                                                                                                                                      >  2:1140950/1‑150 (MQ=255)
                 tATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGTAACTGCCAGGGCtttttt                                                                                                                                    <  1:936909/150‑5 (MQ=255)
                 tATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGAtt                                                                                                                                    >  2:744549/1‑150 (MQ=255)
                 tATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGAtt                                                                                                                                    >  2:731082/1‑150 (MQ=255)
                 tATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGCTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGAtt                                                                                                                                    <  2:296555/150‑1 (MQ=255)
                  aTGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAGAGCCCTGGTAACTGCCAGGGCtttttta                                                                                                                                   <  1:1404415/150‑6 (MQ=255)
                    gTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGTAACTGCCAGGGCttttttatt                                                                                                                                 >  2:591449/1‑143 (MQ=255)
                    gTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGGAGTTAGCAGGGGTTGATTACt                                                                                                                                 >  1:1344793/1‑150 (MQ=255)
                    gTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACt                                                                                                                                 <  2:565093/150‑1 (MQ=255)
                     ttttGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACtt                                                                                                                                >  2:141234/1‑150 (MQ=255)
                      tttGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACttt                                                                                                                               <  1:290492/150‑1 (MQ=255)
                          tCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGc                                                                                                                           <  1:1038938/150‑1 (MQ=255)
                           ctgcctgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCt                                                                                                                          >  1:527955/1‑150 (MQ=255)
                            tgcctgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTa                                                                                                                         >  2:916852/1‑150 (MQ=255)
                             gcctgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTaa                                                                                                                        >  1:286185/1‑150 (MQ=255)
                                  ctctGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAGGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAattatt                                                                                                                   >  1:781231/1‑150 (MQ=255)
                                   tctGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATACAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAAttatta                                                                                                                  >  2:1324658/1‑150 (MQ=255)
                                    ctGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTAc                                                                                                                 <  1:287148/150‑1 (MQ=255)
                                    ctGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTAc                                                                                                                 >  2:963231/1‑150 (MQ=255)
                                      gggATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCCGGGCTTGATTACCTTGAGCGAAttaataatg                                                                                                               >  2:373187/1‑147 (MQ=255)
                                      gggATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTc                                                                                                               <  1:1126187/150‑1 (MQ=255)
                                         ttCGCTGGGGCGGGAATTTTTTTGCCTATTTTGCATTGTTGGTTAGAAAGTATGCAATTCTATTAATTTAAATATGTTGATTAAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAAc                                                                                                            <  2:958665/149‑1 (MQ=255)
                                          tCGCTGGGGCGGGCATTTTTTTTCCTATTTTGGATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACa                                                                                                           <  2:638828/150‑1 (MQ=255)
                                          tCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACa                                                                                                           <  1:1140950/150‑1 (MQ=255)
                                                 ggcgggcATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCg                                                                                                    <  1:522013/150‑1 (MQ=255)
                                                       cATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGAcgcg                                                                                              >  2:662195/1‑150 (MQ=255)
                                                       cATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGAcgcg                                                                                              >  2:1121159/1‑150 (MQ=255)
                                                       cATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGAcgcg                                                                                              >  1:1213991/1‑150 (MQ=255)
                                                                gCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCgg                                                                                     >  1:627709/1‑150 (MQ=255)
                                                                  cTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCa                                                                                   <  2:348388/150‑1 (MQ=255)
                                                                    aTTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCTGACGCGCAGTAGCGGCaga                                                                                 >  2:547575/1‑150 (MQ=255)
                                                                    aTTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCaga                                                                                 >  1:696790/1‑150 (MQ=255)
                                                                       ttGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAg                                                                                                                                                    >  2:233538/1‑80 (MQ=255)
                                                                       ttGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAg                                                                                                                                                    <  1:233538/80‑1 (MQ=255)
                                                                         gCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCat                                                                            >  1:414742/1‑150 (MQ=255)
                                                                          cATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatg                                                                           <  2:1095023/150‑1 (MQ=255)
                                                                              gttgGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGAtgt                                                                       >  1:1282556/1‑150 (MQ=255)
                                                                               ttgGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatgatg                                                                        >  1:1237484/1‑148 (MQ=255)
                                                                               ttgGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatgatg                                                                        <  2:1237484/148‑1 (MQ=255)
                                                                               ttgGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGAtgtg                                                                      >  1:250620/1‑150 (MQ=255)
                                                                                 ggTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGtt                                                                    >  2:6087/1‑150 (MQ=255)
                                                                                  gTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTg                                                                   <  1:65469/150‑1 (MQ=255)
                                                                                     aGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCg                                                                <  1:731082/150‑1 (MQ=255)
                                                                                       cAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTa                                                              <  2:1380370/150‑1 (MQ=255)
                                                                                       cAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTa                                                              >  2:666055/1‑150 (MQ=255)
                                                                                         aGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATg                                                            <  2:286185/150‑1 (MQ=255)
                                                                                               ccATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTTATTATTAGTTAACTGGTTGAGGTAGCGCATGAGGGTAATAGTCTTTATGTGTTGACGTATGacctct                                                      <  2:791240/150‑4 (MQ=255)
                                                                                               ccATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGaccacc                                                      <  2:1213991/150‑1 (MQ=255)
                                                                                                  ttCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGc                                                   <  2:964982/150‑1 (MQ=255)
                                                                                                    cGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCgg                                                 >  1:918568/1‑150 (MQ=255)
                                                                                                    cGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCgg                                                 >  1:911137/1‑150 (MQ=255)
                                                                                                    cGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCgg                                                 >  2:555720/1‑150 (MQ=255)
                                                                                                     gATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGc                                                <  2:918568/150‑1 (MQ=255)
                                                                                                     gATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGc                                                >  2:580768/1‑150 (MQ=255)
                                                                                                     gATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGc                                                <  2:1049051/150‑1 (MQ=255)
                                                                                                        gAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCAcc                                             <  1:873889/150‑1 (MQ=255)
                                                                                                        gAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTAATCAACAGGTTGCGGACGCGCAGGAGCGGAAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCAcc                                             >  1:1367880/1‑150 (MQ=255)
                                                                                                         aaTTTTAATATGTTTATTCAAAGATGAAATAAAAAAGCACTGGCAGTTACAAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGAGCAGGAGCGGCAGAGGAATGATGTGTTGCCGTATGAACACCTGCGGCACCt                                            >  2:79509/1‑150 (MQ=255)
                                                                                                         aaTTATAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTAACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCt                                            <  2:1344793/150‑1 (MQ=255)
                                                                                                              tAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTAcc                                       <  1:1397420/150‑1 (MQ=255)
                                                                                                              tAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTAcc                                       <  1:744549/150‑1 (MQ=255)
                                                                                                                 tatGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTc                                    <  1:507554/150‑1 (MQ=255)
                                                                                                                  atGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCg                                   <  2:627709/150‑1 (MQ=255)
                                                                                                                     ttGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAg                                <  1:1253397/150‑1 (MQ=255)
                                                                                                                     ttGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAg                                <  1:466499/150‑1 (MQ=255)
                                                                                                                            aaaGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAg                         <  1:1121159/150‑1 (MQ=255)
                                                                                                                            aaaGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAg                         >  1:1058082/1‑150 (MQ=255)
                                                                                                                                   aaaTAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGgcg                  >  1:620168/1‑150 (MQ=255)
                                                                                                                                      tAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGa                                  <  1:1304350/131‑1 (MQ=255)
                                                                                                                                      tAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGa                                  >  2:1304350/1‑131 (MQ=255)
                                                                                                                                      tAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTg               >  2:390534/1‑150 (MQ=255)
                                                                                                                                       aatcaaGCCCTGGTAACTGCCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCg                                                                <  2:1213553/96‑1 (MQ=38)
                                                                                                                                       aatcaaGCCCTGGTAACTGCCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCg                                                                >  1:1213553/5‑100 (MQ=38)
                                                                                                                                        atcaaGCCCTGGTAACTGCCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGcc             >  2:1381004/4‑150 (MQ=255)
                                                                                                                                        aaaaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGcc             <  1:615672/150‑1 (MQ=255)
                                                                                                                                        aaaaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGcc             <  2:475334/150‑1 (MQ=255)
                                                                                                                                        aaaaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGcc             <  2:1115420/150‑1 (MQ=255)
                                                                                                                                        aaaaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGcc             >  2:1245240/1‑150 (MQ=255)
                                                                                                                                            aGCCCTGGTAACTGCCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTc         >  2:42767/1‑150 (MQ=255)
                                                                                                                                            aGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTc         >  2:591198/1‑150 (MQ=255)
                                                                                                                                               ccTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACCCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACt      <  2:781231/150‑1 (MQ=255)
                                                                                                                                                  gtcAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACtgtg   >  1:1272147/3‑150 (MQ=255)
                                                                                                                                                  ggCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACtgtg   <  1:598781/150‑1 (MQ=255)
                                                                                                                                                  ggCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACtgtg   <  1:6087/150‑1 (MQ=255)
                                                                                                                                                   gCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGa  <  1:555720/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATATCATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGA  >  NC_000913/1141000‑1141296

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: ‑