Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,858 G→C 6.8% intergenic (‑100/+15) yecJ ← / ← azuC DUF2766 family protein/acid‑inducible small membrane‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,8580GC 6.8% 249.9 / 8.5 74intergenic (‑100/+15)yecJ/azuCDUF2766 family protein/acid‑inducible small membrane‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base C (2/3);  total (37/37)
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.08e-01

CAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCG  >  NC_000913/1987710‑1988004
                                                                                                                                                    |                                                                                                                                                  
caaCCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCg                                                                                                                                                    <  2:401226/149‑1 (MQ=255)
    cacaTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGc                                                                                                                                                >  1:182837/1‑149 (MQ=255)
     acaTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCa                                                                                                                                               <  4:143614/149‑1 (MQ=255)
        tCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGt                                                                                                                                            <  4:525975/149‑1 (MQ=255)
        tCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGt                                                                                                                                            >  3:201060/1‑149 (MQ=255)
         cAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGtt                                                                                                                                           <  2:35963/149‑1 (MQ=255)
            aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGttttt                                                                                                                                        >  4:482201/1‑149 (MQ=255)
             aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGtttttt                                                                                                                                       <  2:31080/149‑1 (MQ=255)
               ccgaTTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACAC‑AAAAAACCTGCTCCGGCAGGtttttttgt                                                                                                                                    <  2:203400/146‑4 (MQ=255)
               ccAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACAC‑AAAAAACCTGCTCCGGCAGGtttttttgt                                                                                                                                    <  4:191806/149‑4 (MQ=255)
                cAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGtt                                                                                                                                    <  1:268308/149‑1 (MQ=255)
                cAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGtt                                                                                                                                    >  3:192678/1‑149 (MQ=255)
                cAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGtt                                                                                                                                    >  3:233701/1‑149 (MQ=255)
                         aTCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACa                                                                                                                           >  2:31491/1‑149 (MQ=255)
                             gCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATAtt                                                                                                                       <  1:31491/149‑1 (MQ=255)
                               aTTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGc                                                                                                                     <  4:334444/149‑1 (MQ=255)
                                 tCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCt                                                                                                                   <  3:143940/149‑1 (MQ=255)
                                  cAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTg                                                                                                                  <  1:218760/149‑1 (MQ=255)
                                    gCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGc                                                                                                                <  3:28108/149‑1 (MQ=255)
                                    gCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGc                                                                                                                >  4:497468/1‑149 (MQ=255)
                                          gggACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACg                                                                                                          >  4:547301/1‑149 (MQ=255)
                                            gACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTc                                                                                                        >  4:159910/1‑149 (MQ=255)
                                                  gTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                           >  3:267692/1‑140 (MQ=255)
                                                  gTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                           <  4:267692/140‑1 (MQ=255)
                                                        tCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAg                                                                                                                                              <  1:550266/99‑1 (MQ=255)
                                                        tCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAg                                                                                                                                              >  2:550266/1‑99 (MQ=255)
                                                          cTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                          >  1:573040/1‑149 (MQ=255)
                                                          cTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                          >  2:174063/1‑149 (MQ=255)
                                                              tAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                                      >  2:593123/1‑149 (MQ=255)
                                                              tAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                                      <  1:593123/149‑1 (MQ=255)
                                                                gCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGtttttt                                                                                                                                       >  4:133106/1‑98 (MQ=255)
                                                                gCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGtttttt                                                                                                                                       <  3:133106/98‑1 (MQ=255)
                                                                gCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGt                                                                                    >  1:322669/1‑149 (MQ=255)
                                                                  aTTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTa                                                                                  >  2:127568/1‑149 (MQ=255)
                                                                   ttttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                                      <  1:143679/144‑1 (MQ=255)
                                                                   ttttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                                      >  2:143679/1‑144 (MQ=255)
                                                                   ttttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTAt                                                                                 <  1:71895/149‑1 (MQ=255)
                                                                    tttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGGttttt                                                                                                                                        >  1:90128/1‑93 (MQ=38)
                                                                    tttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGGttttt                                                                                                                                        <  2:90128/93‑1 (MQ=38)
                                                                    tttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTAtt                                                                                <  3:159910/149‑1 (MQ=255)
                                                                        aaaGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGttttt                                                                                                                                        >  4:100036/1‑89 (MQ=255)
                                                                        aaaGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGttttt                                                                                                                                        <  3:100036/89‑1 (MQ=255)
                                                                                    gtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAgg                                                                >  1:556321/1‑149 (MQ=255)
                                                                                     tagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGaa                                                                                                  >  1:376694/1‑114 (MQ=255)
                                                                                     tagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGaa                                                                                                  <  2:376694/114‑1 (MQ=255)
                                                                                     tagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGa                                                               >  3:494194/1‑149 (MQ=255)
                                                                                       gCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAg                                                                 <  1:298631/145‑1 (MQ=255)
                                                                                       gCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAg                                                                 >  2:298631/1‑145 (MQ=255)
                                                                                               acaTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAg                                                     >  4:101954/1‑149 (MQ=255)
                                                                                                 aTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGtt                                                   >  4:165852/1‑149 (MQ=255)
                                                                                                           aTAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAAc                                                                                                                            >  4:525153/1‑66 (MQ=255)
                                                                                                           aTAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAAc                                                                                                                            <  3:525153/66‑1 (MQ=255)
                                                                                                               tCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtttt                                                           <  1:130302/127‑1 (MQ=255)
                                                                                                               tCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtttt                                                           >  2:130302/1‑127 (MQ=255)
                                                                                                               tCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGtt                                                   <  2:322162/135‑1 (MQ=255)
                                                                                                               tCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGtt                                                   >  1:322162/1‑135 (MQ=255)
                                                                                                                cacaTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTcc                                    >  4:213994/1‑149 (MQ=255)
                                                                                                                      cACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAt                              <  4:201060/149‑1 (MQ=255)
                                                                                                                        ccGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAttt                            <  3:165852/149‑1 (MQ=255)
                                                                                                                           tcatcaGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGtt                         <  1:127568/149‑1 (MQ=255)
                                                                                                                           tcatcaGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGtt                         <  2:556321/149‑1 (MQ=255)
                                                                                                                            catcaGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTta                        <  2:573040/149‑1 (MQ=255)
                                                                                                                                aGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACa                    >  3:211917/1‑149 (MQ=255)
                                                                                                                                aGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACa                    >  2:529083/1‑149 (MQ=255)
                                                                                                                                    cacaAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCAc                                       >  2:314136/1‑126 (MQ=255)
                                                                                                                                    cacaAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCAc                                       <  1:314136/126‑1 (MQ=255)
                                                                                                                                       aaaaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTa             >  3:267455/1‑149 (MQ=255)
                                                                                                                                        aaaaaaCCTGCTCCGGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAt            >  1:59920/1‑149 (MQ=255)
                                                                                                                                         aaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAtt           <  3:316018/149‑1 (MQ=255)
                                                                                                                                         aaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAtt           <  4:233701/149‑1 (MQ=255)
                                                                                                                                         aaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAtt           <  4:192678/149‑1 (MQ=255)
                                                                                                                                            aaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTc        >  4:441822/1‑149 (MQ=255)
                                                                                                                                            aaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTc        <  3:497468/149‑1 (MQ=255)
                                                                                                                                              ccTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCtt      <  4:267455/149‑1 (MQ=255)
                                                                                                                                                  ccGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTgcg  >  1:150195/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCG  >  NC_000913/1987710‑1988004

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