Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,988,106 (C)6→8 intergenic (‑147/‑116) azuC ← / → yecR uncharacterized protein AzuC/lipoprotein YecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,988,1001.C100.0% 221.3 / NA 60intergenic (‑141/‑122)azuC/yecRuncharacterized protein AzuC/lipoprotein YecR
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (36/24);  total (36/24)
*NC_0009131,988,1002.C88.3% 191.9 / 18.6 60intergenic (‑141/‑122)azuC/yecRuncharacterized protein AzuC/lipoprotein YecR
Reads supporting (aligned to +/- strand):  ref base . (4/3);  new base C (32/21);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

GCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAG‑‑‑CCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACTCTTCTGGCTGGACCATG  >  NC_000913/1987964‑1988224
                                                                                                                                         ||                                                                                                                              
gCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAG‑‑‑ccc                                                                                                                            >  1:8686/1‑140 (MQ=255)
gCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCCTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAG‑‑‑ccc                                                                                                                            >  1:190301/1‑140 (MQ=255)
  aCTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAG‑‑‑ccccc                                                                                                                          <  1:957940/140‑1 (MQ=255)
    tCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAG‑‑‑ccccccc                                                                                                                        <  2:741309/140‑2 (MQ=255)
        aCCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAATAAACAG‑‑‑cccccccctat                                                                                                                    >  2:934331/1‑135 (MQ=255)
             ttGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTg                                                                                                                 <  1:203150/140‑1 (MQ=255)
             ttGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTg                                                                                                                 <  1:934331/140‑1 (MQ=255)
              tGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTCTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGa                                                                                                                <  1:809599/140‑1 (MQ=255)
              tGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTCTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGa                                                                                                                >  2:809599/1‑140 (MQ=255)
                       ggTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGa                                                                                                                >  1:810532/1‑131 (MQ=255)
                       ggTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGa                                                                                                                <  2:810532/131‑1 (MQ=255)
                       ggTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAg                                                                                                       >  1:768529/1‑140 (MQ=255)
                        gtgtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGa                                                                                                      >  1:689727/1‑140 (MQ=255)
                        gtgtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGa                                                                                                      >  2:834422/1‑140 (MQ=255)
                          gtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATg                                                                                                    >  2:376967/1‑140 (MQ=255)
                                  cTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGaaa                                                                                            >  2:42242/1‑140 (MQ=255)
                                   tttGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAAt                                                                                           >  1:426531/1‑140 (MQ=255)
                                                    aCAACCATCACACAAATCAATCTTTATGTGCTACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACACATGGCAAGAAATAGCGCCTGCCAGGCGTc                                                                          >  2:203311/1‑140 (MQ=255)
                                                      aaCCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCtt                                                                        >  2:611979/1‑140 (MQ=255)
                                                      aaCCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCtt                                                                        >  2:8213/1‑140 (MQ=255)
                                                       aCCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCttt                                                                       >  1:49448/1‑140 (MQ=255)
                                                       aCCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCttt                                                                       >  1:894492/1‑140 (MQ=255)
                                                        ccATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCtttt                                                                      >  2:689082/1‑140 (MQ=255)
                                                        ccATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCtttt                                                                      >  1:894839/1‑140 (MQ=255)
                                                           tCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCg                                                                   >  1:991292/1‑140 (MQ=255)
                                                            cacaAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGC‑‑CCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGc                                                                 >  1:637756/1‑140 (MQ=255)
                                                                     caatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCg                                                        <  2:470383/140‑1 (MQ=255)
                                                                     caatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCg                                                        <  1:42242/140‑1 (MQ=255)
                                                                     caatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCg                                                        <  1:419230/140‑1 (MQ=255)
                                                                     caatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCg                                                        <  2:894839/140‑1 (MQ=255)
                                                                     caatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCg                                                        <  2:894492/140‑1 (MQ=255)
                                                                      aatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgc                                                       >  2:222184/1‑140 (MQ=255)
                                                                       atcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgca                                                      >  1:203186/1‑140 (MQ=255)
                                                                       atcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgca                                                      <  1:689082/140‑1 (MQ=255)
                                                                             aTGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCAt                                                             >  2:84016/1‑127 (MQ=255)
                                                                             aTGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCAt                                                             <  1:84016/127‑1 (MQ=255)
                                                                               gtgATACAAATCACATAAATTCCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTGCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTCTCGCAGCACTGTa                                              <  2:188898/140‑1 (MQ=255)
                                                                               gtgATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTa                                              <  2:49448/140‑1 (MQ=255)
                                                                                  aTACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTa                                              <  2:826319/137‑1 (MQ=255)
                                                                                  aTACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTa                                              >  1:826319/1‑137 (MQ=255)
                                                                                    aCAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAAcgcg                                         <  2:8686/140‑1 (MQ=255)
                                                                                        aTCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTaaa                                     >  1:359947/1‑140 (MQ=255)
                                                                                          cacaTAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGC‑‑CCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAAc                                            <  2:77916/130‑1 (MQ=255)
                                                                                          cacaTAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGC‑‑CCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCACCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAAc                                            >  1:77916/1‑130 (MQ=255)
                                                                                             aTAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGt                                >  2:974134/1‑140 (MQ=255)
                                                                                               aaaTACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGC‑‑CCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCt                             <  2:637756/140‑1 (MQ=255)
                                                                                                 aTACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCt                             <  1:511329/140‑1 (MQ=255)
                                                                                                  tACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCttt                           >  2:545786/1‑140 (MQ=255)
                                                                                                  tACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCttt                           >  1:40186/1‑140 (MQ=255)
                                                                                                       cTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCtt                      >  1:7709/1‑140 (MQ=255)
                                                                                                       cTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCtt                      >  1:287524/1‑140 (MQ=255)
                                                                                                        tttAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTa                     <  2:203186/140‑1 (MQ=255)
                                                                                                        tttAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCCACAGACGGCAAGAAATAGCGACTGCCAGGCGTCTTTTCCGGCCATTTTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTa                     >  2:943699/1‑140 (MQ=255)
                                                                                                          tAATGTTATAAAAATGATAATCAAAAAACAGC‑‑CCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTActc                  >  1:21391/1‑140 (MQ=255)
                                                                                                           aaTGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTActc                  >  2:814499/1‑140 (MQ=255)
                                                                                                             tGTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACtctt                <  1:853041/140‑1 (MQ=255)
                                                                                                              gTTATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACtcttc               >  2:114300/1‑140 (MQ=255)
                                                                                                               ttATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACtctt                <  1:326744/138‑1 (MQ=255)
                                                                                                               ttATAAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACtctt                >  2:326744/1‑138 (MQ=255)
                                                                                                                 ataAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACTCTTctgg            >  2:732600/1‑140 (MQ=255)
                                                                                                                  taAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTTGTGGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACTCTTctgg            >  2:1003844/1‑140 (MQ=255)
                                                                                                                  taAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTa                     <  1:638592/130‑1 (MQ=255)
                                                                                                                  taAAAATGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTa                     >  2:638592/1‑130 (MQ=255)
                                                                                                                   aaaaaTGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACTCTTctggct          <  1:545786/140‑1 (MQ=255)
                                                                                                                    aaaaTGATAATCAAAAAACAGCCCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACTCTTctggctg         <  1:611979/140‑1 (MQ=255)
                                                                                                                         gATAATCAAAAAACAGC‑‑CCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACTCTTCTGGCTGGACCAt   >  1:317005/1‑140 (MQ=255)
                                                                                                                          aTAATCAAAAAACAGC‑‑CCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACTCTTCTGGCTGGACCATg  <  2:21391/140‑1 (MQ=255)
                                                                                                                                         ||                                                                                                                              
GCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAG‑‑‑CCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTTCTCTTACTCTTCTGGCTGGACCATG  >  NC_000913/1987964‑1988224

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