Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,490,177 C→T 5.5% A2T (GCA→ACA)  argD ← bifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,490,1770CT 5.5% 241.3 / 3.0 73A2T (GCA→ACA) argDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base T (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

CCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGG  >  NC_000913/3490036‑3490325
                                                                                                                                             |                                                                                                                                                    
ccGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCa                                                                                                                                               >  2:35213/1‑149 (MQ=255)
 cgcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCAt                                                                                                                                              <  4:31265/149‑1 (MQ=255)
    gAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGAt                                                                                                                                           <  2:445551/149‑1 (MQ=255)
      aaTCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCa                                                                                                                                         >  2:250074/1‑149 (MQ=255)
       aTCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCAc                                                                                                                                        <  1:466995/149‑1 (MQ=255)
              tACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTa                                                                                                                                 >  4:407253/1‑149 (MQ=255)
              tACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTa                                                                                                                                 <  2:448394/149‑1 (MQ=255)
               aCTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTAc                                                                                                                                >  2:296991/1‑149 (MQ=255)
                    ttgccatgcTGATCCCAGATTCGGCTGACCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGAGTAATTGCTGTTTGTTTAATTGCCATCTCATGATCACCCTGTTACGAATa                                                                                                                           >  2:482004/1‑149 (MQ=255)
                      gccttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATaaa                                                                                                                         <  2:583321/149‑1 (MQ=255)
                      gccttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATaaa                                                                                                                         >  2:363220/1‑149 (MQ=255)
                             tGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAAtgt                                                                                                                  <  1:198791/149‑1 (MQ=255)
                             tGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAAtgt                                                                                                                  >  4:77946/1‑149 (MQ=255)
                             tGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAAtgt                                                                                                                  >  2:43293/1‑149 (MQ=255)
                              gATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAAtgtg                                                                                                                 >  3:279784/1‑149 (MQ=255)
                                tCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGaa                                                                                                               >  2:520872/1‑149 (MQ=255)
                                 cccAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACaaa                                                                                                                     >  1:557409/1‑142 (MQ=255)
                                 cccAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACaaa                                                                                                                     <  2:557409/142‑1 (MQ=255)
                                 cccAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACaaa                                                                                                                     >  3:229035/1‑142 (MQ=255)
                                 cccAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACaaa                                                                                                                     <  4:229035/142‑1 (MQ=255)
                                  ccAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAAt                                                                                                             >  3:345063/1‑149 (MQ=255)
                                     gATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATtat                                                                                                          <  2:382234/149‑1 (MQ=255)
                                       ttCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATaa                                                                                                        <  1:425036/149‑1 (MQ=255)
                                        tCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAAc                                                                                                       <  1:520872/149‑1 (MQ=255)
                                         cGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAAcc                                                                                                      <  3:65247/149‑1 (MQ=255)
                                               cccTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACaaaa                                                                                                <  1:35213/149‑1 (MQ=255)
                                                ccTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAat                                                                                               >  1:552987/1‑149 (MQ=255)
                                                  ggACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAatat                                                                                             >  2:479146/2‑149 (MQ=255)
                                                   gACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATg                                                                                            >  3:160990/1‑149 (MQ=255)
                                                       ttttACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGTCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATa                                                                                        >  1:145113/1‑149 (MQ=255)
                                                       ttttACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGTCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATa                                                                                        >  3:464065/1‑149 (MQ=255)
                                                       ttttACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATa                                                                                        >  3:189628/1‑149 (MQ=255)
                                                               gAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCa                                                                                <  3:407253/149‑1 (MQ=255)
                                                               gAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCa                                                                                >  4:294206/1‑149 (MQ=255)
                                                               gAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCa                                                                                <  1:250074/149‑1 (MQ=255)
                                                                aaTAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCAc                                                                               >  2:576687/1‑149 (MQ=255)
                                                                     aCTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGACGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAAt                                                                          <  4:177495/149‑1 (MQ=255)
                                                                      ctctGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATg                                                                         >  3:370368/1‑149 (MQ=255)
                                                                       tctGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATgg                                                                        <  1:488648/149‑1 (MQ=255)
                                                                       tctGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATgg                                                                        <  2:471300/149‑1 (MQ=255)
                                                                           ccGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAAt                                                                    <  4:79556/149‑1 (MQ=255)
                                                                             ggAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCa                                                                  >  1:292121/1‑149 (MQ=255)
                                                                                       tCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCg                                                        >  4:136229/1‑149 (MQ=255)
                                                                                          gCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATg                                                     >  1:48878/1‑149 (MQ=255)
                                                                                            aGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCa                                                   >  2:471689/1‑149 (MQ=255)
                                                                                             ggATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAg                                                  <  4:21535/149‑1 (MQ=255)
                                                                                             ggATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAg                                                  <  1:576687/149‑1 (MQ=255)
                                                                                                tCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGaa                                               <  1:463021/149‑1 (MQ=255)
                                                                                                tCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGaa                                               >  3:355099/1‑149 (MQ=255)
                                                                                                 cACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGaaa                                              <  2:522122/149‑1 (MQ=255)
                                                                                                 cACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGaaa                                              <  1:436055/149‑1 (MQ=255)
                                                                                                 cACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGaaa                                              >  4:352725/1‑149 (MQ=255)
                                                                                                    ttCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTa                                           >  4:444415/1‑149 (MQ=255)
                                                                                                       aTCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGcc                                        <  1:295384/149‑1 (MQ=255)
                                                                                                       aTCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGcc                                        >  3:47939/1‑149 (MQ=255)
                                                                                                        ttgAAAGTCGCGCGGGTAATTGCAGTTTGTTCAATTGCCATCTCATCATCACCCTGTTACGCATATACAAATGTGAAATTATAACCACAAAATATGAATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCca                                       >  2:192978/3‑149 (MQ=255)
                                                                                                        tCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCca                                       >  3:336559/1‑149 (MQ=255)
                                                                                                           aaaGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCcagca                                    <  3:294206/149‑1 (MQ=255)
                                                                                                           aaaGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCcagca                                    <  2:48878/149‑1 (MQ=255)
                                                                                                              gTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGtt                                 <  2:253323/149‑1 (MQ=255)
                                                                                                              gTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGtt                                 <  2:292121/149‑1 (MQ=255)
                                                                                                                   gtgtgtAATTGCTGTTTGTTCAATTGCCATCACATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGt                            <  2:261180/147‑1 (MQ=255)
                                                                                                                    cgTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTc                           <  1:574521/149‑1 (MQ=255)
                                                                                                                          aTTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGtt                     <  3:77946/149‑1 (MQ=255)
                                                                                                                          aTTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGtt                     <  1:43293/149‑1 (MQ=255)
                                                                                                                                gtttgttCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTaa               <  2:101558/149‑1 (MQ=255)
                                                                                                                                gtttgttCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTaa               >  2:151915/1‑149 (MQ=255)
                                                                                                                                      tCAATTGTCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATAc         <  2:145113/149‑1 (MQ=255)
                                                                                                                                      tCAATTGTCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATAc         <  4:464065/149‑1 (MQ=255)
                                                                                                                                        aaTTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTa                                           <  2:103205/113‑1 (MQ=255)
                                                                                                                                        aaTTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTa                                           >  1:103205/1‑113 (MQ=255)
                                                                                                                                          ttGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTc     <  4:370368/149‑1 (MQ=255)
                                                                                                                                          ttGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTc     <  3:111801/149‑1 (MQ=255)
                                                                                                                                             ccATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCtgg  >  1:23580/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGG  >  NC_000913/3490036‑3490325

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