Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 985,855 A→T 7.3% intergenic (‑146/+39) aspC ← / ← ompF aspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913985,8550AT 7.3% 254.5 / 8.1 96intergenic (‑146/+39)aspC/ompFaspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (53/36);  new base T (3/4);  total (56/40)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

ATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTT  >  NC_000913/985708‑986001
                                                                                                                                                   |                                                                                                                                                  
aTGACGATGTTACATTCTGGTTACAGTAGTGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAATTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCaaa                                                                                                                                                   <  2:361323/149‑1 (MQ=255)
aTGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCaaa                                                                                                                                                   >  2:1207968/1‑149 (MQ=255)
  gACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGt                                                                                                                                                 >  1:1078355/1‑149 (MQ=255)
     gAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACAAAAGTCCt                                                                                                                                              >  1:900423/1‑149 (MQ=255)
       ggTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGt                                                                                                                                            >  1:832314/1‑149 (MQ=255)
        gTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGCCCAAAGTCCTGtt                                                                                                                                           >  2:311614/1‑149 (MQ=255)
          tCCATTATGGTTACAGAAGTGAAGTCAGAAATCAGGGTAACGGGAGATTTACAAAAATCCAACTATTACTTATGAAAACGAATGCTGTTTTTTCAATCCGTGAGAATGCTCTTGAAGGTATAAAAAAAACAGGACCAACGTCCTGtttt                                                                                                                                         <  1:1292960/149‑1 (MQ=255)
          tCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGtttt                                                                                                                                         <  1:72216/149‑1 (MQ=255)
          tCCATTATGGTTACAGAAGGGAAGTCCGCAATCAGGGTACCGGGAGATATACAAACTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACTTTGGTCCTGtttt                                                                                                                                         <  2:410925/149‑1 (MQ=255)
           ccATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACTTTGGTCCTGttttt                                                                                                                                        <  2:497108/149‑1 (MQ=255)
            cATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGtttttt                                                                                                                                       >  1:870859/1‑149 (MQ=255)
             aTTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACTTTGGTCCTGttttttt                                                                                                                                      <  2:1506884/149‑2 (MQ=255)
                 tGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAAGTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCGCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGct                                                                                                                                  >  2:1384409/1‑148 (MQ=255)
                 tGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCa                                                                                                                                  <  2:55131/149‑1 (MQ=255)
                   gTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAtt                                                                                                                                <  1:1274334/149‑1 (MQ=255)
                   gTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAtt                                                                                                                                >  1:1094593/1‑149 (MQ=255)
                    ttACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACGATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGAAttt                                                                                                                               >  2:644417/1‑149 (MQ=255)
                       cAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAAc                                                                                                                            <  1:105411/149‑1 (MQ=255)
                       cAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAAc                                                                                                                            <  2:379022/149‑1 (MQ=255)
                       cAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACa                                                                                                                           <  1:737100/149‑1 (MQ=255)
                        aGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACa                                                                                                                           >  1:804773/1‑149 (MQ=255)
                            gggAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAaga                                                                                                                       >  1:762079/1‑149 (MQ=255)
                             ggAAGTCCGCTATCAGGGTAACGTGAGATTTACAAAATTCCAACTATTAATGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAagag                                                                                                                      <  1:133770/149‑1 (MQ=255)
                             ggAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACACAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCTGCATTTAACAAagag                                                                                                                      >  2:721756/1‑149 (MQ=255)
                                 gTCCGCTATCAGGGTAACGGGAGAATTACAAAATTCCAACTATTACTGATGAAACCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtgt                                                                                                                  <  2:1493087/149‑1 (MQ=255)
                                        aTCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTAGGCAATTAACAAAGAGGTGTGGTAttt                                                                                                           >  2:720812/1‑148 (MQ=255)
                                          cAGGGTCACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGTTGTGCTATTAGa                                                                                                         >  2:296082/1‑149 (MQ=255)
                                            gggTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAAc                                                                                                       >  2:1253496/1‑149 (MQ=255)
                                              gTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTg                                                                                                     <  1:429662/149‑1 (MQ=255)
                                                  cGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCTGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaa                                                                                                 >  2:300135/1‑149 (MQ=255)
                                                   gggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCTGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaaa                                                                                                <  1:998755/149‑1 (MQ=255)
                                                   gggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaaa                                                                                                <  1:650672/149‑1 (MQ=255)
                                                     cagaTTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCATGTTATTTCGGCATTTAACAAAGAGATCTGCTATTAGAACTGGTAAACg                                                                                              >  2:700911/2‑149 (MQ=255)
                                                         tttACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAcc                                                                                         <  1:1369365/149‑1 (MQ=255)
                                                           tACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCTTTTAACAAAGAGGTGTGCTATTAGAACTGGTTAACGATAccc                                                                                        >  2:124895/1‑149 (MQ=255)
                                                            aCAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaaa                                                                                                <  1:134950/140‑1 (MQ=255)
                                                            aCAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaaa                                                                                                >  2:134950/1‑140 (MQ=255)
                                                            aCAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACAAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCca                                                                                       <  2:832314/149‑1 (MQ=255)
                                                              aaaaTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcaca                                                                                     >  1:1216576/1‑149 (MQ=255)
                                                              aaaaTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTTCTATTAGAACTGGTAAACGATACCCACAg                                                                                    >  2:1313456/1‑149 (MQ=255)
                                                               aaaTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAg                                                                                    >  1:754398/1‑149 (MQ=255)
                                                                     cAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGGTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACgg                                                                              >  2:1547085/1‑149 (MQ=255)
                                                                     cAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAActg                                                                              >  2:412036/1‑147 (MQ=255)
                                                                      aaCTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGgt                                                                             <  2:1078355/149‑1 (MQ=255)
                                                                       aCTATTACTGATGTAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAAAGgtg                                                                            <  2:903630/149‑1 (MQ=255)
                                                                        cTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGACCTGTTTTTTCGGCATTTAACAATGAGGTTTGCTATTATAACTGGAAAACGATACCAACAGCAACGgagt                                                                           >  1:245670/1‑149 (MQ=255)
                                                                         tATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTgtc                                                                          <  1:300135/149‑1 (MQ=255)
                                                                         tATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCCGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGgagta                                                                          >  1:521797/1‑148 (MQ=255)
                                                                           tttcTGATGAAAACGCACGCTGTTTTTGAAAGAAGTGAGATTGCTCTGGAAGGTATAAAAAATACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTgtcgt                                                                        <  2:1082165/146‑1 (MQ=255)
                                                                            tACTGATGAAAACGCAGGCTGTTTTTGCAACACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGACCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTgtcgta                                                                       >  2:1253865/1‑148 (MQ=255)
                                                                                gatgaACACGCCTGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAATGTATAAAAAAAACAGGACCAAAGTCCTGTTTATTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGaa                                                                   <  2:792351/149‑1 (MQ=255)
                                                                                 atgaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCAAAATCAACGGAGTCGTCTGAAc                                                                  >  1:1460779/1‑149 (MQ=255)
                                                                                       aCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGc                                                            >  1:738685/1‑149 (MQ=255)
                                                                                         gCAGGCTGTTTTTGCAAGACGTGAGATTGCACTCGAAGGTATAAAAAAAACAGGACCAAAGAACTGTTTTTTCGGCATATAACAAAGAGCTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGCTGTCGTCTGAACATACTccc                                                          >  2:1191581/1‑149 (MQ=255)
                                                                                           aGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAg                                                        >  1:553528/1‑149 (MQ=255)
                                                                                           aGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAg                                                        >  2:1312426/1‑149 (MQ=255)
                                                                                             gCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAACACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACATAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACCGCAACGGTGGCGTCTGAACCTACGCCCAGtt                                                      <  2:1112414/149‑1 (MQ=255)
                                                                                             gCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGtt                                                      >  2:592616/1‑149 (MQ=255)
                                                                                                 tttttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAATTCATGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttgtt                                                  <  2:1040220/149‑1 (MQ=255)
                                                                                                    tggCAAGACGTGAGATTGCTCTGGAAGGTTTAAAAAAAACAGGACCAAAGGCCTGTTTTTTCGGCATTGACCAAAGAGGTGCGCTATTCTACCTGGTCAACGATACCCACAGCACCGGTGTCGTCTGAACCTACGCCCAGTTTGATGTc                                               <  2:1010227/147‑1 (MQ=255)
                                                                                                     ttcatgACGTTAGTTTGCGATGGATGGTTTAAAAAATACAGGTGCATTGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTTGAACTGGTAAACGATACCCACTGCTACGGTGTCGTCTGAAGGTAGGCCCAGTTTTTTGTCa                                              <  2:116813/144‑1 (MQ=255)
                                                                                                      gCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAg                                             >  2:31313/1‑149 (MQ=255)
                                                                                                       cAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGa                                            >  1:757364/1‑149 (MQ=255)
                                                                                                        aaGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAATCGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGat                                           >  2:1545792/1‑148 (MQ=255)
                                                                                                        aaGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGaa                                           <  2:553528/149‑1 (MQ=255)
                                                                                                            cGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGa                                       >  2:1330166/1‑149 (MQ=255)
                                                                                                             gTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGTTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGAt                                      <  1:721756/149‑1 (MQ=255)
                                                                                                             ctGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGAt                                      <  2:1434574/148‑1 (MQ=255)
                                                                                                               gagaTTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAATAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCt                                    >  2:1534968/1‑149 (MQ=255)
                                                                                                               gagaTTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCt                                    >  1:1204801/1‑149 (MQ=255)
                                                                                                               gagaTTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCt                                    <  1:788689/149‑1 (MQ=255)
                                                                                                                 gaTTGCTCTGGAAGTTATAAAAAAAACAGGCACAAAGTACTGTTTTTTCGGCATTTAACAAAGAGGTGTGATATTATAACGGGTAAACTATACCCACAGAAACGGTGTCGTCTGAACCTAATCCCAGTTTGTTGTCAGAATCGATCTgg                                  <  2:900423/149‑1 (MQ=255)
                                                                                                                   tttctctGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGtt                                <  1:777334/146‑1 (MQ=255)
                                                                                                                      ctctGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCCCAGCAACGGTGTCGTCGGTACCTACGCCCAGTTTGTTGTCAGAATCGAGCTGGTtgag                             >  1:1010426/1‑148 (MQ=255)
                                                                                                                       tctGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGGTTTTTCGGCATATAACAAAGAGGTGTGCTATTAGAACTGGTAATCGATACCCCCAGCAACGGTGTCGTCTGACCCTACGCCCAGTTTGTTGTCAGAATCGCTCTGGTtgatg                            >  1:474392/1‑149 (MQ=255)
                                                                                                                        gtGGAAGGTATAAAAAAAACAGTACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGATATTAGAACTGGTAAACGAGACCCACAGCACCGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatta                           <  1:508574/148‑3 (MQ=255)
                                                                                                                         tGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTTGTAgattat                          >  2:338994/1‑149 (MQ=255)
                                                                                                                           gAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGt                        >  1:450318/1‑149 (MQ=255)
                                                                                                                            aaGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTa                       <  2:870859/149‑1 (MQ=255)
                                                                                                                               gTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGTCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCACAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGta                    >  2:1486045/1‑148 (MQ=255)
                                                                                                                                 ataAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAGCCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCaa                  >  1:142272/1‑149 (MQ=255)
                                                                                                                                 ataAAAAAAACAGGACCAAAGTCCTGGTTTTTCGGCATTTAACAAAGAGGTGTCCTATTAGAACTGGTAAACTATACCCACAGCAACGGTGTCGTCTGAACCTACGCCGAGTTTGTTGACAGATTCGATCTGGTTGATGATGTAGTaaa                  >  1:1099127/1‑149 (MQ=255)
                                                                                                                                  taAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAAc                 <  1:296082/149‑1 (MQ=255)
                                                                                                                                  taAAAAAAACAGGACCAAAGTCCGGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTGTGAACCAAAGCCCAGTTTGGTGTCAGAATAGATCTGGTTGATGATGTAGTCAAc                 >  2:329762/1‑149 (MQ=255)
                                                                                                                                   cgaaaaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACCAAGAGGTGTGCTATTAGAACTGGGAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGCTCTGGTTGATGATGTAGTCAACa                >  1:901278/3‑149 (MQ=255)
                                                                                                                                   aaaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAAAAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACa                >  1:496693/1‑149 (MQ=255)
                                                                                                                                    gaaaaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCACCGGTGTCGTCTGATCCTACGCACAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTcatcat               >  1:1194636/2‑149 (MQ=255)
                                                                                                                                    aaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACAt               >  1:507982/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATa              >  2:1405369/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGGTGATGATGTAGTCAACATa              >  2:1133664/1‑149 (MQ=255)
                                                                                                                                      aaaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACACACAGCAACGGTGTCGTCAGAACCTACGCCCAGTTTGTTGTCAGAAGCGATCTGGTTGATGATGTAGTCAACATAg             >  2:471677/1‑149 (MQ=255)
                                                                                                                                      aaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAg             <  2:762079/149‑1 (MQ=255)
                                                                                                                                          aCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTgg         <  2:757364/149‑1 (MQ=255)
                                                                                                                                          aCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTgg         >  1:605469/1‑149 (MQ=255)
                                                                                                                                          aCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTgg         <  1:587513/149‑1 (MQ=255)
                                                                                                                                              gACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACAt     >  1:464064/1‑149 (MQ=255)
                                                                                                                                              gACCAAAGGCCTGTTTTTTCGGCATTTACCAAAGAGGTTTGCTATTAGACCTGGTAAACGATACCCACAGCAACGGCGTCGGCTGAACCTACGCCCAGTTTGTTGTCAGTATCGATCTGGATGATGATGTAGTCAACATAGGTGGACAt     <  1:575717/149‑1 (MQ=255)
                                                                                                                                                 cAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGATTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGAGGACATGtt  >  1:1252733/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTT  >  NC_000913/985708‑986001

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: