Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,288,518 A→G N252S (AAC→AGC)  nrfA → cytochrome c552 nitrite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,288,5180AG97.5% 269.1 / ‑2.7 81N252S (AAC→AGC) nrfAcytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (36/42);  total (38/43)
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.49e-01

CGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGCAAACT  >  NC_000913/4288390‑4288653
                                                                                                                                |                                                                                                                                       
cGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCt                                                                                                                                       >  1:222985/1‑131 (MQ=255)
     aGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACg                                                                                                                         >  2:130896/1‑140 (MQ=255)
      gTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGc                                                                                                                        <  1:1191534/140‑1 (MQ=255)
        gCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCa                                                                                                                      <  1:878420/140‑1 (MQ=255)
               ggAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGa                                                                                                               <  2:1069495/140‑1 (MQ=255)
                  gTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAg                                                                                                            >  2:1172085/1‑140 (MQ=255)
                   tATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAgc                                                                                                           >  1:1251685/1‑140 (MQ=255)
                    aTTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAgcg                                                                                                          <  1:1330590/140‑1 (MQ=255)
                    aTTACTTCGACGGCAAAAACAAAGCGGTTAAATTCACGTCGGATGACGCCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAgcg                                                                                                          <  2:467156/140‑1 (MQ=255)
                     ttACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAgcgc                                                                                                         <  2:229555/140‑1 (MQ=255)
                        cTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGGCTACCTCCCTGTCGAAAACGCAAATGCTGAAAGCgcagc                                                                                                      >  2:499384/1‑140 (MQ=255)
                        cTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCgcagc                                                                                                      >  2:864860/1‑140 (MQ=255)
                        cTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCgca                                                                                                        >  2:771032/1‑138 (MQ=255)
                        cTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCgca                                                                                                        <  1:771032/138‑1 (MQ=255)
                        cTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACGAGCTCCCGGTCGAATACGCCAATGCTGAAAGCgcagc                                                                                                      >  1:503474/1‑140 (MQ=255)
                         ttCGACGGCAAAAACAAAGCGGCTAAATTCCCGTGCGATTACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCgcagca                                                                                                     <  2:552771/140‑1 (MQ=255)
                           cGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATTACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCAcc                                                                                                   <  1:605713/140‑1 (MQ=255)
                           cGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACACAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGcccc                                                                                                   <  2:889975/140‑1 (MQ=255)
                           cGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCAcc                                                                                                   <  2:1169113/140‑1 (MQ=255)
                           cGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCAcc                                                                                                   >  2:1130871/1‑140 (MQ=255)
                           cGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCAcc                                                                                                   >  2:274027/1‑140 (MQ=255)
                           cGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCAcc                                                                                                   >  2:904778/1‑140 (MQ=255)
                           cGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCAcc                                                                                                   >  2:1526642/1‑140 (MQ=255)
                            gACGGCAAAAACAAACCGGTTAGATTCCCTTGGGATGTCGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCAccc                                                                                                  <  1:701110/140‑1 (MQ=255)
                             aCGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCAcc                                                                                                   <  2:17288/138‑1 (MQ=255)
                             aCGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCAcc                                                                                                   >  1:17288/1‑138 (MQ=255)
                             aCGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCg                                                                                                 >  2:769313/1‑140 (MQ=255)
                             aCGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCg                                                                                                 <  1:904778/140‑1 (MQ=255)
                             aCGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCg                                                                                                 >  1:944349/1‑140 (MQ=255)
                                 cAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAat                                                                                             >  1:1040454/1‑140 (MQ=255)
                                  aaaaaCAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGGCGAAAACACCAATGTTGAAAGAGCAGCACCCGGAata                                                                                            >  1:552771/1‑140 (MQ=255)
                                   aaaaCAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAatat                                                                                           <  2:873386/140‑1 (MQ=255)
                                   aaaaCAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAatat                                                                                           <  1:178111/140‑1 (MQ=255)
                                   aaaaCAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAatat                                                                                           <  2:802484/140‑1 (MQ=255)
                                   aaaaCAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAatat                                                                                           <  1:130896/140‑1 (MQ=255)
                                                       ccGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCAt                                                                       >  1:757821/1‑140 (MQ=255)
                                                         gTGGGATGACGTCATGAAATGCGGAATTATGGAGCAGTATTACGACAACATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTc                                                                     >  2:1132131/1‑140 (MQ=255)
                                                         gTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTc                                                                     <  2:262823/140‑1 (MQ=255)
                                                         gTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTc                                                                     <  1:1130871/140‑1 (MQ=255)
                                                         gTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTc                                                                     >  1:1178342/1‑140 (MQ=255)
                                                             gATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACgg                                                                 <  1:861458/140‑1 (MQ=255)
                                                              aTGACTGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGt                                                                <  2:802125/140‑1 (MQ=255)
                                                              aTGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGt                                                                <  2:374654/140‑1 (MQ=255)
                                                              aTGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGt                                                                >  2:85178/1‑140 (MQ=255)
                                                                      aTGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAaacaa                                                        >  2:1537937/1‑140 (MQ=255)
                                                                       tGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAaacaac                                                       <  2:1178342/140‑1 (MQ=255)
                                                                       tGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCACTGACTGGACTACCTCCCTGTCGTAAACGCCGAAGCCGAAACGGCCGCATCCGCAATATGAAATCTGGGCAGCGGGCATTCACGGTAAAatcaac                                                       <  2:376954/140‑1 (MQ=255)
                                                                        gAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACg                                                      >  2:1426765/1‑140 (MQ=255)
                                                                         aaaGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGt                                                     <  2:264936/140‑1 (MQ=255)
                                                                           aGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGa                                                   <  1:788087/140‑1 (MQ=255)
                                                                              cGAAAATATGGAGCAGTATTACGACAAAATTTCCTACTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCt                                                <  2:503474/140‑1 (MQ=255)
                                                                               gAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAGTATGAAACCTGGACAGCGGGCATTCACGGTaaaa                                                            >  2:1331775/1‑127 (MQ=255)
                                                                               gAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAGTATGAAACCTGGACAGCGGGCATTCACGGTaaaa                                                            <  1:1331775/127‑1 (MQ=255)
                                                                                aaaaTATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGt                                              <  1:1172085/140‑1 (MQ=255)
                                                                                    tatGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCg                                          <  1:1537937/140‑1 (MQ=255)
                                                                                                   cGACAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCaa                           >  1:1372638/1‑140 (MQ=255)
                                                                                                     aCAAAATTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGAACTATATCGAGTGCCATATGCCaaaa                         >  1:629204/1‑140 (MQ=255)
                                                                                                        aaaTTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGTACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTg                      >  2:261495/1‑140 (MQ=255)
                                                                                                        aaaTTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAACAACAACGCGACCTGTAGCGACAGCCATATGCCAAAAGTg                      >  2:652269/1‑140 (MQ=255)
                                                                                                        aaaTTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTg                      >  1:1441939/1‑140 (MQ=255)
                                                                                                        aaaTTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTg                      >  2:215177/1‑140 (MQ=255)
                                                                                                        aaaTTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTg                      >  1:143314/1‑140 (MQ=255)
                                                                                                        aaaTTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTg                      >  2:1415485/1‑140 (MQ=255)
                                                                                                        aaaTTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTg                      >  1:786733/1‑140 (MQ=255)
                                                                                                          aTTGCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCa                    <  1:1135959/140‑1 (MQ=255)
                                                                                                             gCCTTCTCTGACTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGaa                 >  2:1462816/1‑140 (MQ=255)
                                                                                                                       aCTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGTACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGc       <  1:261495/140‑1 (MQ=255)
                                                                                                                       aCTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGc       <  1:1462816/140‑1 (MQ=255)
                                                                                                                       aCTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGc       <  2:1287809/140‑1 (MQ=255)
                                                                                                                       aCTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGc       <  2:944349/140‑1 (MQ=255)
                                                                                                                       aCTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGc       <  1:1415485/140‑1 (MQ=255)
                                                                                                                       aCTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGc       <  1:864860/140‑1 (MQ=255)
                                                                                                                       aCTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGc       <  2:1040454/140‑1 (MQ=255)
                                                                                                                       aCTGGACTAGCTCCCTGTCGAAAACCCCACTGCTGAACGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGc       <  1:499384/140‑1 (MQ=255)
                                                                                                                        cTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGCa      <  2:1372638/140‑1 (MQ=255)
                                                                                                                        cTGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACACCGAAGGCa      <  2:757821/140‑1 (MQ=255)
                                                                                                                        cTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGCa      <  1:85178/140‑1 (MQ=255)
                                                                                                                         tGGACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGCaa     >  2:1100277/1‑140 (MQ=255)
                                                                                                                           gACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGc       <  1:1420670/136‑1 (MQ=255)
                                                                                                                           gACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGc       >  2:1420670/1‑136 (MQ=255)
                                                                                                                           gACTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGCAAAc   >  1:239886/1‑140 (MQ=255)
                                                                                                                            aCTAGCTCCCTGTCGAAAACGCCTATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGCAAAct  >  2:310104/1‑140 (MQ=255)
                                                                                                                            aCTAGCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGCAAAct  >  2:1249038/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                       
CGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAGTGCAGAACGCCGAAGGCAAACT  >  NC_000913/4288390‑4288653

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